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PREFACE 



This report is one in a 'series on the University of Miami's thermal 
pollution models. Much of the background, formulations, solutions tech- 
niques and applications of these models has been summarized in a three- 
volume report by Lee and Sengupta (1978). These mathematical models 
were developed by the Thermal Pollution Croup at the University of Miami 
and were funded by the National Aeronautics and Space Administration, 
The primary aim was to have a package of mathematical models which has 
general application in predicting the thermal distribution of once-through, 
power plant heated discharge to the aquatic ecosystem. The joint effort 
was planned so that the calibration and verification of these models de- 
pend on simultaneous remote sensing and ground-truth data acquisition 
support. The concept is the development of adequately calibrated and 
verified models for direct prediction of thermal dispersion by the user 
communities. The intended user communities include the utility companies 
and the regulatory agencies at the federal and state level. 

The purpose of the present effort is to further verify these models 
at widely different sites using minimal calibration, and then to provide 
the program code and user's manual to the Environmental Protection 
Agency (EPA) for future users. 

Two sites chosen were Anclote Anchorage on the west coast of 
Florida and Lake Keowee in South Carolina. The free-surface model was 
applied to Anclote Anchorage and the rigid-lid to Lake Keowee. Two 
data acquisition trips, one in the summer and the other in the winter, 
were carried out at each site. The acquisition was a collaborated effort 
jointly by the UM, EPA, NASA-KSC (Kennedy Space Center) and the 
corresponding utility company personnels. 

The two-year project consisted of two phases. During Phase I, 
the individual model was modified and calibrated to fit the corresponding 
site; then the model was verified against the remote sensing isotherms 
and in-situ measured velocity and temperature data under both summer 
and winter conditions. In general, the computed isotherms were compar- 
able with isotherms based on remote sen<tnq. In Phase II, source pro- 
grams of both models were documented and transferred to the EPA, and 
the user's guides were prepared to familiarize the potential users. The 
results of the two-year effort are summarized in a set of three final re- 
ports. 



ABSTRACT 



To assess the environmental impact of waste heat dlsDOsal by power 
plant operation into natural water bodies, mathematical models are essen- 
tial, especially for predictive studies. The Thermal Pollution Croup at 
the University of Miami has developed a package of models for this pur- 
pose. A joint effort with EPA, NASA, Duke Power Company and Florida 
Power Company was conducted to verify these models with remote sensed 
IR data and in-situ measurements. 

The free-surface model, presented in this volume, is for tidal 
estuaries and coastal regions where ambient tidal forces play an important 
role in the dispersal of heated water. The model is time dependent, three 
dimensional and can handle irregular bottom topography. The vertical 
stretching coordinate is adopted for better treatment of kinematic condition 
at the water surface. The results include surface elevation, velocity and 
temperature* 

The model has been verified at the Anclote Anchorage site of Flo* :da 
Power Company. Two data bases at four tidal stages for winter and summer 
conditions were used to verify the model. Differences between measured and 
predicted temperatures are on an average of less than 1°C. 
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SECTION 1 
INTRODUCTION 



BACKGROUND 

The problem of disposing of the waste heat produced as a result of 
generation of electrical energy, whether by fossil-fuel or by nuclear-fuel, 
is a dominant consideration in making power production compatible with 
ecological concerns. For every unit of energy converted to electricity, 
approximately two units are rejected in the form of waste heat. The ul- 
timate heat s>nk for the heat removed by the condenser-cooling-water 
system is V . arth's atmosphere. The cooling water taken from natural 
or man-made water bodies is circulated once through the condenser, and 
the heated water is discharged back to the same bodies, which usually 
are lakes, » : vers, estuaries, or coastal waters* Eventually, the heat is 
transferred to the atmosphere through evaporation, radiation and conduc- 
tion over relatively large areas at the air-water interface. 

The use of natural water bodies as an intermediate means for dispo- 
sal of waste heat must take into account the effect upon the environment 
of the circulation and temperature rises produced in the receiving water. 
The rate of oxyge^ consumption by aquatic species increases with rising 
water temperature, however, the ability of water to hold dissolved oxygen 
decreases with nsing water temperature. There are possibilities of im- 
pairment of biological functions of fish and of breaking important links 
in the food chain. The lethal effects of thermal pollution are sometimes 
obvious; the sublethal effects on hydrobiological systems and waste assi- 
milation capacities are not easy to foresee unless interactive hydrothermal, 
chemical and biological studies are conducted in an integrated fashion. 

Accurate understanding of hydrothermal behavior of the receiving 
water bodies is an important factor in a power plant system for the fol- 
lowing reasons: 

1. To provide a priori information about the nature and extent of ther- 

mal impact on the aquatic life forms. 

2. To analyse the circulation pattern of the receiving water body so 
that recirculation between intake and outlet and consequent decrease 
in cooling efficiency can be minimized. 

3. To assess the thermal impact on the aquatic life forms existing in the 
receiving ecosystem so that post-operational remedies can be done to 



reduce the hazards. 

It is therefore apparent that not only environmental but planning 
and designing interests also are at stake* 

The above-mentioned objectives can only be met by having large 
data sets over the entire discharge flowfleld under various extreme hydro- 
logical and meteorological conditions. Measurement for temperature and 
velocity made over the affected domain could be used to develop maps for 
velocity and temperature distributions. These in-situ measurements can 
serve for diagnostic and monitoring purposes under limited circumstances; 
however, they are not relevant for predictive objective. The physical 
modeling is useful for hydrodynamical behavior studies when properly 
ve/ified; it is generally expensive and time-consuming. Frequently, 
physical dimensions necessitate distortion in the model, making exact 
dynamic similitude impossible. Mathematical modeling is tractable and 
predictive; it has become a powerful means in simulating complex environ- 
mental flows. 

In order to establish the rationale of model formulation, the physical 
mechanisms underlying the heat dispersion from a heated discharge need 
to be outlined. The following mechanisms govern the heat dfspersal. 

1. Entrainment of ambient fluid into the thermal discharge. 

2. Buoyant spreading of discharged heated effluent. 

3. Diffusion by ambient turbulence. 
4* Interaction with ambient currents. 

5. Heat loss to the atmosphera through air-water interface. 

The first four mechanism? redistribute heat and momentum in the re- 
ceiving water body. The last mechanism eventually transfers heat to the 
atmosphere. It has been customary, therefore, to make assumptions and 
approximations which enable the model solvable. For example, the ambient 
turbulence is considered by assuming the eddy viscosity and diffusivity 
dependent on the mean velocity field. There assumptions necessitate 
careful calibration of models to assure reliability. Large data bases are 
needed for proper calibration, especially for three-dimensional, time-de- 
pendent models. Remote sensing is the only available method of obtaining 
large synoptic data bases. Sengupta et al. (1975) has discussed the need 
for remotely sensed data for adequate development of time-dependent, 
hydrothermal models. 

UNIVERSITY OF MIAMI MODELS 

Critical evaluation of mathematical models used for thermal pollution 
analysis has been made by Dunn et al. (1975). They compared the per- 



formance of various models in predicting a standard data base. A gene- 
ral conclusion that can be made from their analysis is that, though some 
models may perform well under certain conditions, a generalized three- 
dimensional model which accounts for wind, current, tide, bottom topo- 
graphy and diverse meteorological conditions is yet to be developed. 
One of the first three-dimensional models was by Waldrop and Farmer 
(1973, 1974 a, b) . This model was essentially a free-surface formula- 
tion. One of the first three-dimensional models which adequately accounts 
for bottom topography and comprehensive meteorological conditions was a 
rigid-Hd model developed by Sengupta and Lick (1974 a, 1976). They 
used a vertical stretching to convert the variable depth, a concept cus- 
tomarily adopted by numerical weather forecasting peoples. 

The thermal pollution research team at the University of .Miami has 
for the past several years been developing a package of three-dimensional 
mathematical models which could have general application to problems of 
power plant heated discharge to the aquatic ecosystem. The primary 
motivation behind the effort was to develop a series of models with mini- 
mal restrictive assumptions, enabling applications to diverse basin confi- 
gurations and to various driving forces of ambient flow. The effort is 
closely integrated with simultaneous remote sensing and ground-truth 
data acquisition support. Our aim is to develop adequately calibrated 
and verified models for production purpose; that is, for direct applica- 
tion by the user communities. The user communities are the power in- 
dustries and the regulatory agencies like the Environmental Protection 
Agency and the Nuclear Regulatory Commission. 

For the time-dependent, free-surface model which is the main concern 
of this volume, the governing equations are the incompressible Navier- 
Stokes equations, conservation of mass, energy and an equation of state. 
The conventional hydrostatic and Boussinesq approximations are made. 
The mean velocity field closure is assumed in defining eddy viscosity 
coefficient. The boundary conditions at water-land interfaces are no- 
slip, no-normal velocity and adiabatic conditions. At the air-water inter- 
face, wine! stress and heat transfer coefficients are specified. At open- 
to-sea boundaries, conditions are specified for surface elevation and tem- 
perature. Otherwise, normal derivatives of temperature and velocity are 
equal to zero. Initial condition is assumed to be equilibrium; that is 
"cool start." Conditions at intake and out*** are completely specified in 
space and time. 

The features of the UM's free-surface model can be summarized as 
follows: 

1. Gravity, earth's rotation, nonlinear inertia and bottom friction terms 
are included in the hydromechanic part. 

2. The driving forces include wind, tide, river outflow and power plant 
intake and discharge. 

3. Convection, diffusion and heat transport at water surface are 

3 



included in the thermodynamic part. 

4. The a system of coordinate, devised by a meteorologist for avoiding 
difficulties at the free surface, is incorporated, 

5. It predicts spatial and temporal variation of surface elevation due 
to tide surge and wind set-up. 

6. It has means for graphically representing velocity and temperature 
fields. 

This free-surface model is general enough to be applied with 
minor modifications to a large variety of sites, such as lakes, rivers, 
estuaries, tidal inlets and coastlines. This model has been successfully 
applied to Lake Okeechobee for hydrodynamic study, Biscayne Bay for 
dispersion study and Hutchinson Island for hydrothermal study. This 
model can simulate the receiving water body in the far-field and de- 
tailed features of thermal plumes and mixing zones in the near-field. 

DESCRIPTION OF ANCLOTE ANCHORAGE 

The UM's thermal pollution team has been studying the Anclote 
Anchorage on the west central coast of Florida near the town of Tarpon 
Springs (Figure I) since the summer of 1978. The Anclote power plant, 
operated by the Florida Power Corporation, has two 515 MW, oil-fired 
electrical generating units* The once-through cooling water for the 
two units is to be drawn from the Anclote River by six pumps deliver- 
ing a total of 930,000 gpm. After a temperature rise of 6. 1°C, this 
water is diluted with a flow of 1,060,000 gpm at ambient temperature to 
reduce the temperature rise to 2. 8°C at the outlet to Anclote Anchorage. 
This mixing is done in a 1250 meter-long man-made channel leading to 
an outlet. 

The Anclote Anchorage (Figure 2) consists of shallow coastal water 
separated from the Culf of Mexico by a series of barrier island parallel 
to the coast line. The Anchorage has a relatively unrestricted exchange 
of water with the Culf through natural channels to the north and to the 
south of the Anclote Keys. Depth ranges from 0. 3 to 3. 6 m, with a 
mean of 1.8 m. Shallow regions of less than 0.6 m comprise approximately 
5 km in length and 6 km in average width. Currents in the Anchorage 
are tidal- and wind-driven, with the tide entering from the south stronger 
than that from the north. 

Prior to power plant construction, the Marine Science Institute of 
the University of South Florida (USF) was contracted by the Florida 
Power Corporation (FPC) to investigate the possible environmental im- 
pact of the plant operation. The USF's Anclote Environmental Project 
beginning in 1970 was comprehensive in nature. Much of their efforts 
was to obtain a detailed picture of the Anclote environment prior to 
alternation of the environment by power plant construction and operation. 



Valuable background data were collected in a variety of areas , including 
general physical characteristics of the Anclote River and Anchorage, 
suspended sediments, turbidity, temperature, salinity, water quality, 
seagrasses, benthic invertebrates and fishes. All information was ap- 
plied to preconstruction plant design in order to minimize environmental 
impact. This application resulted in major changes in intake and dis- 
charge channel dredging, outfall design and thermal dishcarge character- 
istics. The same team also monitored the undergoing changes of the 
environment during and after the dredging operations for the intake and 
discharge channels. 

Unit I of the Anclote Power Plant commenced operation during the 
fall of 1974, after which a post-operation ecological monitoring program 
for evaluating the nature and degree of thermal impact was carried out 
and maintained by the USF. At the time of this study, the planned 
Unit 2 was still pending permission. 

Since the nature of thermal impact at Anclote Anchorage is not 
clear, it is decided that a joint effort by UM, NASA-KSC, EPA and FPC 
to use a three-dimensional model with support of remotely sensed data 
and in-situ data for calibration and verification may be appropriate to 
study this thermal impact at Anclote. It is for this purpose "That the 
existing free-surface model was adapted to the Anclote site, as we have 
done in its application to Biscayne Bay. As the original program was 
developed for well-mixed shallow coastal waters, except for some modi- 
fication to accommodate the tidally influenced Anclote River flow, the 
program was ready for calibration. We shall account for this application 
in a later section. 



SECTION 2 
CONCLUSIONS 



A numerical simulation of the hydrothermal characteristics of a well- 
mixed, shallow, coastal water body at Anclote Anchorage under the effect 
of waste heat disposal by a power plant in the summer and winter situa- 
tions is presented* The model takes into consideration the effects of 
Coriolis force, wind, tide, bottom topography, power plant intake and 
discharge, river outflow and surface heat transfer. Results obtained 
with the model have been verified with in-situ measurements and IR data. 
Reasonable agreement has been obtained* From the experiences and re- 
sults of these simulation runs, the conclusions may be summarized as 
follows: 

1. Inputs for the description of open boundaries, discharge, wind and 
heat transfer must be found experimentally from field data for ac- 
curate hydrothermal predictions. 

2. The shape of the thermal plume was dominated by the stage of the 
tidal cycles, as clearly exhibited in the plots of IR and calculated 
isotherms. 

3. Tide plays a main role in the Anchorage as a driving force; its 'in- 
fluence from the south is stronger than that from the north. A 
dividing line is observed. On the northern side of this line, the 
water flows north, while on the south side it flows south. At this 
ridge line, the water has minimal transverse motion. The location 

of this line varies with time as the tide from both ends is not in phase. 

4. It is important to impose correct boundary conditions, especially the 
correct tidal functions on the south and north boundaries, for ob- 
taining a good prediction of the thermal plume. 

5. Wind does appear to be an important external force affecting surface 
currents. It is found that the surface currents are close to the 
wind direction when wind speed is in excess of 15 mph. 

6. An estimation of the Rossby Number revealed that the nonlinear 
inertia terms can be safely neglected for Anclote site. This is be- 
cause the ratio of inertia force to the Coriolis force is small compared 
with unity. 

7. The recirculation of cooling water should be prevented by proper 
design of intake and discharge location, since the heated water 



recirculating back to the intake will cause a reduction in the effi- 
ciency of the power plant. 

Numerically the model behaves very well for both summer and winter 
simulation runs. The model is able to include very shallow water depths. 
A problem associated with the analysis is computing the response of the 
thermal plume to very strong winds such as hurricane, since the formula 
used in this study to estimate the vertical eddy viscosity coefficient does 
not involve wind speed* A brief investigation was made regarding the 
effects of hurricane-force winds. It was found that currents become un- 
realistically large in such cases unless the vertical eddy viscosity coeffi- 
cient was increased with wind speed. This general problem of computing 
the response of shallow coastal water to very strong winds requires addi- 
tional research. 



SECTION 3 
RECOMMENDATIONS 



One of the difficulties encountered in model verification with the 
remote sensing temperature field was the lack of measured current data 
of comparable accuracy. We believe that, at the verification state, the 
accuracy of velocity calculation needs to be assessed in order to clear 
ambiguities about the limitations of comparison of calculated surface iso- 
thermal maps with the IR imagery mosaics. 

In the initial phase of mixing / the plume shape is governed by the 
volume of discharge, the geometry of outlet and the initial temperature 
difference between the discharge and receiving writers. Consequently, 
for predicting the details of the near-field thermal plume a finer grid is 
needed near the outlet. Thus, a combination of grid structures would be 
desired for the computation with a course grid for far-field and a fine 
grid for the near-field. This, however, would result in higher computa- 
tional costs. 

For any non-reactive, dissolved, chemical constituent, the governing 
transport equation is completely analogous to the temperature equation 
of the present model. Thus, without much endeavor, the model can be 
extended to include the dispersal of these constituents. 

For situations of deep and stratified water, the bouyancy effects 
are important factors in shaping the plume. Therefore, a coupled system 
of momentum and energy equations should be considered as the basis of 
formulation. 



SECTION 4 
MATHEMATICAL FORMULATION AND MODEL DESCRIPTION 

GENERAL BACKGROUND ON FREE-SURFACE (TIDAL) MODELS 

Recent concern about the ecological future of our estuaries and 
coastal waters has generated a need for practical and reliable methods of 
predicting the environmental impact of the widespread use of natural sea- 
shores as industrial zones. The tendency of setting up once-through, 
seawater-cooled, fossil- or nuclear-fueled power generation plants is in- 
creasingly strong, due to the demand for electrical power at competitive 
prices and prohibitive costs of cooling towers and man-made cooling lakes* 
Thus, it is of special importance that mathematical models be developed 
for simulating and predicting thermal pollution. The estuaries, tidal 
inlets and coastal waters which serve as the receiving body for waste 
heat are usually of complex configuration and topography. The flow is 
driven by tide, wind, run off and buoyancy force. Thus, a complete 
hydrothermal model for coastal waters must, in addition to solving the 
three-dimensional equations of mass, momentum and energy, include the 
salinity equation to determine closely the local density. A model of this 
type is still far beyond feasibility. Therefore, assumptions and approxi- 
mations must be made to render the system closed and tractable. 

Hindwood and Wallis (1975) have comoil^d a bibliography of 141 
papers concerning computer models for tidal hydrodynamics. In general, 
the output from the hydrodynamic /hydraulic model is recorded on maanetic 
tape. This tape is then entered into the thermal dispersion model. The 
dispersion model operates through successive solution of the finite differ- 
ence equations for the change of temperature (with time and space) due 
to diffusion (eddy mixing), advectlon (velocity transport), and heat 
flux through the water surface. Therefore, it is appropriate to say that 
the hydrodynamic /hydraulic model serves as the backbone to the thermal 
pollution model of the well-mixed water. Commonly, the hydrodynamic 
models are two-dimensional and based on the vertically-integrated equations 
of motion and continuity for an incompressible fluid. Thus, the vertical 
structure of the circulation was not considered. This procedure turned 
out to be sufficient for the investigation of tidal processes. 

Of the two-dimensional models, we should briefly mention the Leen- 
dertse (I967) model and Reid and Bodine (1 969) model for their comprehen- 
siveness and popularity. In Leendertse's model, the unknowns are verti- 
cally-averaged velocities and water level; the bottom friction is in terms 
of Chezy coefficient which is to be calibrated. Reid and Bodine used 



vertically integrated transports and water level as unknowns and Darcy 
friction coefficient. Both models used the space-staggered system which 
gives the simplest scheme consistent with the control volume approach to 
deriving difference equations* However, in the Leendertse model, the 
finite difference equations are expressed in an alternating direction im- 
plicit (ADD form, with two successive time-level operations being executed 
during each time cycle; while in the Reid and Bodine model a time-split 
explicit (leap-frog) method is used for marching forward, the time step 
is limited by numerical stability requirement (cf. Platzman, 1958), less 
than the value niin(Ax, Ay)/ 2gD a ) , in which D mav is the maximum 
depth. max max 

In the context of two-dimensional t»dally-driven flows, we should 
include the finite element method models winch were developed in the last 
few years. This is a method combining finite element for spacial discre- 
tization and finite difference for temporal discretization. The variables 
in discrete element is approximated by simple polynomials whose coeffi- 
cients are expressed in terms of nodal values of the variables and their 
derivatives. Correctly formulated, the physical conservation principles 
are satisfied and, in theory, the element shape and interpolation function 
are quite free as long as certain compactibility conditions are satisfied. 
The freedom of using irregular grid to fit complicated geometry is the 
primary advantage of this method. 

At the present stage, FEM (finite element method) models for tidal 
hydrodynamics are all two-dimensional. Wang's CAFE (Circulation Analy- 
sis by Finite Element) model and Brebbia's shallow water model (in his 
Finite Element Hydrodynamic Problem Orientated Language (FEHPOL) 
package) are both productive. CAFE (Wang, 1978) has a linear triangu- 
lar element for all variables, vertically-integrated transports and water 
level, and has split-time method for time integration. Patridge and 
Brebbia (1976) use a six-node, quadratic triangular element for all vari- 
ables; in this case, they ar^ vertically-averaged velocities and free- 
surface elevation. The 4th order Runge-Kutta (explicit) and Trapezoidal 
Rule (implicit) methods are used for time integration. In general, the 
flexibility of the grid layout, the consistency of FEM formulation and the 
easiness in taking into account the spatially variable properties are the 
advantages. However, the time integration scheme is very much problem- 
oriented, and more fruitful research could lead to a better time integra- 
tion method. It is worthwhile to point out that, in the FDE (finite 
difference method) for two-dimensional hydrodynamics, the explicit me- 
thod had advantage over the implicit method in terms of computer time. 

In the classical approach of two-dimensional computations of tides 
ar,d storm surges in a shallow sea, it is assumed that the velocities are 
uniform over the vertical; however, the vertically uniform velocities are 
of little evidence in the case of those propagation and transport processes 
which are essentially vertically structured. The application of a three- 
dimensional model requires extensive computation on a large number of 
grid points if relatively good resolution is needed. This precludes the 
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use of the implicit method which requires that matrix solutions of a row 
or column of variables be found. In fact, Leendertse et al. (1973) deve- 
loped a three-dimensional model wF ch really uses the explicit leap-frog 
method instead of the implicit methods such as the ADI method Leendertse 
used in his two-dimensional model* 

Leendertse' s three-dimensional model is a vertically-integrated, multi- 
layered model, in which the usual assumptions of hydrostatic pressure, 
imcompressibility and small density variation are made. The interfaces 
are assumed as fixed horizontal planes, while the bottom layer has its 
lower face fit to bottom topography and the top layer has its upper face 
representing the free surface. The top layer has a time-variable height. 
The other layers may be intersected by the bottom and have a height 
which is dependent on the bathymetry. The description of the finite- 
difference equations from the differential equations is accomplished in 
two steps: First, the equations for the layer are derived by vertically 
integrating the variables over the layer thickness, and subsequently, 
finite difference approximations for the layer equation are developed. 
Thus, the vertically-integrated momentum, heat, and salt balance equa- 
tions can be presented for each layer. From the continuity equation, 
the time derivative of the free surface and the vertical velocity at the 
interfaces can be derived. The horizontal pressure gradient- in each 
layer is approximated from hydrostatic equation with layer-averaged 
density, which is tied to salt and heat through the equation of state. 

In Leendertse's model (Leendertse and Liu, 1975), the balance of 
momentum fluxes between the local and convective accelerations, pressure 
gradient, Coriolis force, lateral diffusion, and interfacial stresses is 
accounted for within each layer; likewise, the balance of heat fluxes be- 
tween the time rate, convection, lateral diffusion and cross-layer diffu- 
sion is for each layer; so is the salt balance* The same lateral diffusion 
coefficient is used for both heat and salt. Even postulations are intro- 
duced to express the vertical momentum, salt and heat exchange under 
vertically stable or unstable stratifications; these layer-averaged equa- 
tions clearly point out the foremost problem of his model. There are 
many parameters to be determined and an enormous amount of supporting 
field data is needed. 

Leendertse et al. (1973) discuss the numerical finite difference 
solution scheme in some detail. The explicit leap-frog method is used to 
avoid the difficulties encountered by ADI* The spatial grid structure is 
cell-like with u, v and w components at the center of the corresponding 
normal faces and pressure (p) , density (p), salinity (S) and temperature 
(T) *t the center of the cell. For programming reasons, the bottom must 
be approximated in steps of layer thickness, causing some numerical pro- 
blems at the jumps. But the most troublesome one is the specification of 
boundary conditions for the seaward boundary of the model, since field 
measurements do not provide enough data to describe that boundary 
completely. Numerical experiments involving the simulations of hydrody- 
namic behavior of both Chesapeake Bay and San Francisco Bay are given 
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with a graphical representation of velocity, salinity and temperature results* 

A specific three-dimensional model has been developed by Sundermann 
(1974) and applied to the North Sea. In the model, the usual assumption 
of a hydrostatic pressure and the Boussinesq approximation for the turbu- 
lent Reynolds stresses are made. However, no horizontal turbulent ex- 
change of momentum is retained. To fit the vertical velocity profiles of 
a wind-generated surface current and a compensating countercurrent near 
the bottom, as observed in nature, it was found necessary to include a 
bottom boundary layer that can be modeled by having a vertical eddy 
viscosity dependent on the depth. The spatial discretization is carried 
out by means of a cubic grid net, while the explicit method of Crank- 
Nicolson has been used for approximating the vertical diffusion term. 

The c system of coordinates, devised by N. A. Phillips (1957), was 
used for numerical forecasting. Its advantages are that the kinematic 
boundary conditions at free-surface and bottom are made simple and a 
vertical stretching is used to avoid the difficulty of using regular grid 
net for irregular bathymetry. Briefly, the (x, y, cr, t) system is used 
to replace (x, y, z, t) , and the free surface and bottom are trans- 
formed into a = and or = -1 surfaces. The modified vertical velocity, 
n = da/dt, is zero at both these surfaces, while the actual vertical 
velocity is w = dn/dt at the free surface and w = at the bottom. 

Freeman et al. (1972) introduced ^transformation into their three- 
dimensional free-surface model for wind-generated circulation in a closed 
region. Sengupta and Lick (1 974) incorporated crtransformation in their 
rigid-lid model for wind-driven circulation in lakes* However, the trans- 
formation was used only for the advantage of vertical stretching so that 
a constant vertical grid size, i j, can be used throughout the domain. 
This rigid-lid approximation effectively eliminates surface gravity wave 
by imposing a zero-actual vertical velocity at the water surface* An 
additional equation, a Poisson equation for surface pressure which con- 
tains the rigid-lid condition, can be derived. Therefore, at each time 
step the Poisson equation has to be solved by iterative method. However, 
for large problems, particularly when normal derivative boundary condi- 
tions are used, this can be a time consuming part of the calculation pro- 
cedure. As mentioned, this model is not for tidally-driven coastal flow. 

The present free-surface model is similar to that of Freeman et al. 
in using cr-coordinate system. However, the horizontal turbulent exchange 
of momentum terms is neglected while the vertical turbulent exchange term 
is approximated by using the Dufort-Frankel differencing to ease the con- 
straint on numerical stability due to neglect of the horizontal diffusion 
terms. This model was developed by Carter (1 977) for his study on wind- 
driven flow in Lake Okeechobee. Later, Sengupta et al. (1978) applied 
it to tidal flow in Biscayne Bay. 

GOVERNING EQUATIONS AMD BOUNDARY CONDITIONS , 

Flow of water in estuaries and coastal areas is predominantly hori- 
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zontal. Vertical velocities do occur and are important, as they character- 
ize the vertical circulation. However, the vertical acceleration is extremely 
f^iall in comparison with the gravitational acceleration, so that the vertical 

-celeration is neglected and the vertical equation of motion is then re- 
placed by the hydrostatic assumption* The effect of density variations 
on the inertial and diffusion terms in the horizontal equations of motion 
is neglected. Density variation is retained in the lateral pressure gradi- 
ent terms; that is, the effect of bouyancy is accounted for by allowing 
density variations in the horizontal pressure gradients which influence 
the fluid motion through the horizontal momentum equations. Consequently, 
the continuity equation replaces the mass conservation equation. In the 
estuaries, the density is influenced by the salinity and temperature; in 
this model, the density will be taken dependent only on the temperature. 

For turbulent closure, the Boussinesq approximation for the turbu- 
lent Reynolds stresses is made, *nd lateral dispersion of momentum is 
further assumed unimportant when compared to the rest of the terms in 
the horizontal equations t f motion* The vertical eddy viscosuy, the 
vertical and horizontal e' dy diffusivity of temperature are assumed to be 
constants, although the horizontal eddy diffusivity has orders of magni- 
tude larger than the vertical eddy diffusivity, being due to the much 
larger horizontal characteristic length, L, in comparison with- the vertical 
characteristic length, H. 

With these conditions, we can write the set of governing equations 
expressing the conservation of mass, mementum, and energy in incompres- 
sible flow: 



3u + 3uu + 3uv + 3uw + I 3d . 3ju = 

3t 3x oy 9z o 3x v^z - (1) 



3v L ?uv ^ 5vv ^ 3vw ^ f ^ 1 3p - 3 2 v A 



|E ., 0g = o 



32 vs " (3) 



5H +21 + i^ =0 

3x 3y 3i, (4) 



3T . 3uT . VT . 3wT „ , " : T 



+ 



^V - B v (~4 + r—l) - B ^4 = 



3t j a 3y ' 32 h 3x- ' 3y* v3z- (5] 

o = j(T) (6) 
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where x, y, 2 = Cartesian coordinates positive eastward, northward and 
upward, respectively 

u, v, w = Respective components of velocity in x, y, z direction 

t = Time 

p = Pressure 

p = Density 

T = Temperature 

f = Corioiis parameter 

g = Gravitational acceleration 

A y =* Vertical eddy viscosity 

B. = Horizontal eddy diffusivity 

B = Vertical eddy diffusivity 

The first three equations represent the equations of motion. Equa- 
tion (4) is the equation of continuity, and Equation (5) represents the 
energy equation. The equation of state, Equation (6), expresses the re- 
lation between the density and the temperature. For a more complete 
representation of an estuarine ecosystem, an equation expressing the 
balance of the salts dissolved in the water should be included together 
with balance equations of dissolved substances that are important to the 
analysis* All these equations have the same form as Equation (5)< If 
salinity is included then the density should be related to both salinity 
and temperature. 

Boundary conditions for the above equations must be specified: 

u=v=w=0 at z= -h(x. y) (7) 

\ 



p = 



2H 

3z 



21 

3z 



3t 



x ^v V 

oA„ r 3z = pA„ 



K 



> 



at z = r,(x. v, t) 



iC P B v e 



+ v— - 



3C 



V 



w = 



/ 



(3a) 



(3b) 



(8c) 



(3d) 
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where n = Free-surface elevation above mean water level 

t , t = Surface wind stressei 
x y 

C = Specific heat at constant pressure 

K = Surface heat transfer coefficient 

T = Equilibrium temperature 

T = Surface water temperature 

Also, boundary conditions at closed and open boundaries must be known. 

One major difficulty in the treatment of the free-surface model is 
the fulfillment of kinematic condition at the free surface- The approach 
used in the model is to follow a vertical stretching transformation sug- 
gested by Phillips (1957) and used successively in lake circulation studies 
by Freeman et a!„ (1972). Using this j-transformation, the free surface 
becomes a fixed, flat surface and the variable depth bottom becomes a 
fixed, flat bottom of constant depth* Thus, this method allows easy 
adaptation to various bottom topographies, and a constant vertical grid 
size can be used throughout the domain. 

With respect to the mean water level (MWL) , z is n(x, y, t) at the 
free surface and -h(x, y) at the bottom* The --transformation of the 
vertical coordinate for the free-surface model is obtained by introducing 

z ~ Mx, y, t) _ z - - 

= h(x. y) + n(x. y. t) " H(x. y. t) (9) 

where H = h + r is the total depth measured from local water surface* 
Figure 2 shows the (x, y, j) coordinate system. Note that the value 
of c decreases monotonically from zero at the free surface to minus unity 
at the bottom. The modified vertical velocity, f;, is cieariy zero at both 
free surface and bottom. From Equation (9), the actual vertical velocity, 
w, is related to Q by 

w = (h + n)f. ♦ (j + 1)|~ + s(u^~ * vy r ) (10) 

Differential transformation relationships are required to convert 
equations in (x, y, z) to those in (x, y a). The first derivatives can 
be written as 

( ^ } y.z,t = ( ^ 3 y.a.t ^ H Sx " ,x j 3a (lla) 
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l 3y J x,z,t v 3y'x,a.t H l0 3y 3y'^ (lib) 



riZi -IE. 

l 3z J x,y.t "H 3o (lie) 

where F is any dependent variable and H = h + n is the total depth and 
is independent of a. The continuity equation can be written as 

3x 3y 3a 3t v ; 

By integrating Equation (12) from the bottom to the surface and observing 
the fact that ft vanishes at either face, we obtain the first useful form of 
the continuity equation. 



in - - f r— + 3Hv< 

3t ~ " '-I 3x 3y 



(13) 



The horizontal momentum equations are 

3HU 3Huu 3Hvu + „3£u + (1 + o) |H |l . fvH 

at 3x 3y ?c eta at 

H 3p „„, an . 3ns . A v 3 2 u (14) 



3Hv + 3HUV + 3Hvv + H igv |v 3r, + fuH 

at 3x 3y da aa 3t 

H 32 ,ur 3H , an, . A v 3 2 v 

= " p ay " s H(a a^ + p + H^" (15) 

The energy equation can be written as 



!HT + «ul + «Zl +H ^I + (l + o)g^=B h H(g + |^) + -^0 (16 , 



The coupling of momentum and energy equations may be retained via 
the density which is assumed to be a function of temperature only. An 
empirical formula 
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p = 1.029431 - O.Q0002QT - O.OOOOOST* (17) 

may be used for sea water of salinity 38 parts per thousand. The hydro- 
static equation can be integrated to obtain a diagnostic equation for p. 

p(c) = ~2 H [ a p (a)da ( 18 ) 

Jo 

The continuity equation, Equation (12), can also be integrated from 
the free surface to a to yield the modified vertical velocity a at a plane. 
Thus, a second useful form of Equation (12) yields 

a (a) = u -r-r - u (-^r- + "^rr)da 



1 In 1 f c 

= H 3t"H 

Jo 



3X 3y^ w (19) 



The actual vertical velocity w is recovered through relationship Equation 

(10). In solving parabolic-type partial-differential Equations (14)- (16) 

and differential Equation (13), we choose the initial conditions to be zero 

elevation from MWL and zero velocities. Civen the initial temperature 

field, the boundary conditions are u = v = w = on all solid" surfaces, 

and at a = 0, the wind stresses (t ,t ) are exerted; i.e. 

x y 



A . 


A . 


v 3u 


v 3v 


pH do ~ T x f 


-m r-= t at a 
pn Sa y 



(20a, b) 

The adiabatic conditions are assumed on all solid surfaces, and at a = 0, 
the heat flux is set proportional to the difference between surface tempera- 
ture T and equilibrium temperature T ; i.e. 

pC p B v 3T _ K ,- T i n+ „ - n 

— -£ — — = (T - T ) at a = frM 

H 8c s e s (21) 

The specification of velocities and temperatures at the open boundar- 
ies where the tide enters the basin is more difficult. In the present study, 
the tides outside the basin are given, and the difference of surface eleva- 
tions across the open boundaries is used to determine the normal velocities 
there. The temperature at the open boundaries is set equal to ambient 
temperature. The velocities and temperature at the discharge is specified 
according to plant operation. 

UNCOUPLED SYSTEM 

The foregoing system may be uncoupled so that the hydrodynamic 
mode! for tidal-driven flow is separated from the thermal dispersion model. 
The former system consists of Equations (13), (1U), (15) and (19), while 
the latter contains Equation (16) only. This decoupling amounts to an 
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assumption that the density variations due to the temperature rises 
caused by waste heat dispersion are negligibly small, so that the water 
is of constant density and the energy equation can be uncoupled from 
the continuity and momentum equations. The lateral pressure gradient 
is thus directly related to the gradient of the free surface. 

In the present problem of shallow water flow strongly influenced by 
tide, the thermal dispersion is mainly caused by velocity transport and 
the buoyancy effect is negligible with comparison to the tidal effect. 
Numerically, this decoupling has two obvious advantages. Firstly, the 
time step of the coupled model is the same as that of the hydrodynamic 
model. This time step is limited by the criterion for computational 
stability by the explicit method; namely, At must satisfy the condition 
cAt/As ■< 1, where As represents one of the space intervals in the three- 
dimensional grid and c denotes the maximum characteristic speed. In 
the free-surface model, c is identified with the wave speed and As with 
horizontal spacing; thus, 

At < m * n ( Ax ' A y) 

U - J* H max 

However, the stability analysis of linearized energy equation alone yields 
the thermal time step to be 



^ " [M + M + 2B r_L_ + _L 
1 Ax Ay "Vax 1 Ay 



sH 



where u and v are the maximum particle velocities in x and y direction 
respectively. For the present problem, these time steps are 50 sec and 
400 sec respectively; therefore, in using At =15 sec, we can store the 
hydrodynamic results every 20 steps to matcfi it- = 300 sec used in the 
detached calculation of temperature field. Secondly, the decoupling 
allows one to try for flow solution before it is used for temperature cal- 
culations, and the appropriate flow solution may be used for many tem- 
perature solutions of various initial and boundary conditions of tempera- 
ture. 

The effect of stratification is known to arrest the thermal dispersion. 
For a shallow basin having an insignificant river discharge in comparison 
with tidal flow, which is the case of concern, the basin is nearly well- 
mixed; the density variation is small. Thus, uncoupling is physically 
sound and numerically beneficial. 

COMPUTATIONAL GRID 

To represent the equations in finite difference form, a horizontal 
staggered computing grid system is used. Its plan version, shown in 
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Figure 4, indicates the arrangement of field variables in the x, y plane. 
Figure 5 shows a vertical fluid column subdivided into four layers; each 
has a constant nondimensional thickness, Aa = 0,25. The u, v and w 
velocities are shown at their definition point respectively. Having divided 
the region into cells by a series of grid points which are spaced at dis- 
tances of Ax, Ay, and Ac, the time variable is differenced into increment 
of At such that t * nAt, where n denotes the current time step. 

The water depth, h. ., is specified at full-grid point where both 
horizontal indices, i and l J | are integers; however, since the total depth, 
H = h + n, is needed at half-grid point, where both indices i and j have 
half-integer values, (i + i, j + J}# the specified water depth at full-grid 
points is averaged* The following notations are used to indicate water 
depth at the u-, v- and n-points. 

h i u ,i s * (h i.j + h i.j+i ): h i v .j = * (h i.j + Vi.j' (22a - b} 



^V *% + Vi.j + Vi + Wi> (23 > 



A notation for n at full-grid point is needed for calculating the 
nonlinear inertia terms of the momentum equations; 

E U = * ( Yi + Vi.i + Vi.j-i + Vk' (24) 

where n indicates the present-time surface elevation above the mean 
water level. With Equations (22) to (24) and the calculated water level, 
rv . at present time step, the present total water depths at u-, v- and 
n-points are respectively given by the following expressions. 

U, * =h 4 U ; + i(n n _, , + nf.) 



i.j an i.J + * tn i-l.J + V (25) 



(26) 



H i.j -h i.j + Yj (27) 

In the momentum equations, values of u and v are required at half-grid 
points where values of these variables are not defined and thus, not 
stored. In these cases, values are obtained by linear interpolation be- 
tween the values stored for that variable at the two neighboring points. 
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Therefore, the following average notations are adopted 

n ^ . n 



<j = * (u I+l.j + "ij (28) 

V f.j = Hv i n .i + ^.j+l 3 ( 29) 

"i.j = i(u l.i + "t.f-l' (Ml 



^lj " * ( <J * "5-l.f 1 (33> 

It is clear that notations (28) and (29) give velocities at half-grid point 
while notations (30) and (31) give velocities at full-grid point. No 
superscript n is needed for these notations* These and other variables 
are shown in Figure 6, where u-, v- and n-points are marked. 

FINITE DIFFERENCE EQUATIONS 

In the equations presented here, the space subscripts are i, j, k 
and are all integers, unless otherwise noted. The superscript n refers 
to the time level; n is present time, at which the predictive equations, 
(13) through (16), are used to advance the fields of u, v, n and T to 
new time level n + 1. In this study, first derivatives with respect to 
time are always represented by central time difference, i.e., for any 
variable F, 

F a+1 - f 11 : 1 

l 3t'i.j.k 2At uut J (32) 

where 0(At 2 ) refers to the order of truncation error associated with this 
differencing. The central time scheme is thus of second order accuracy, 
since the neglected terms are terms multiplied by (At 2 ) or a higher order 
of At. This indicates that the finite difference approximation can be 
made to approach the differential by taking sufficiently small At. 

The second derivatives with respect to o in Equations (14) to (16) 
are written in the DuFort-Frankel format, i.e. 

r 3 2 F } n , t.j,k-l i,j, k *i,j,k l,j,k+l 

l 3cr 5 " J i f j.k - (id- (33) 

This expression for second order cr-derivatives when used with the central 
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time difference is known to avoid the inherent instability caused by using 
the common form of second order difference; i.e. 

Hn£)!\ k - F i n i,k-i " g i n ,j,k + F iYk+i 

3a*i,j.k jj^p 

The continuity equation, (13), is the predictive equation for surface 
elevation n, and thus applies at rr point. For convenience, we introduce, 
at level k of (i,j) -column, a variable 

F n = U Hl.iVl,i.k- U i.i u i,i,k , V l,^l v i,Hl,k" V i.i V i.i.k 
k Ax Ay 

where the superscript n on the right-hand-side variables has been dropped. 
Then the finite difference form of the continuity equation can be written as 

n+1 n-1 

Using Simpson's rule, the above equation gives 

n iT = Vj 1 + <*">¥ [F 1 + W 2 + ^3 + ^4 + F 5 ] (34) 

It is to be noted that F_ = since both u and v vanish at the solid bottom. 

Another form of continuity equation, (19), can be used to calculate 
the modified vertical velocity a at each level of the column. That is 

.. . , , n+1 n-1 

_ -(k-l)A a n, . - n, . 

1,J,K H i,i 2At 



+ H £ - {iF l +F 2 + -"- +iF k ] (35) 



where the trapezoidal formula has been used and H. . is from Equation (27); 
i.e. I/J 

H, . =H" =h n . + tP, 
i.J i.] i.] M 
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In the computation, always the present time fl only is needed; therefore, 
no superscript n is required for J2 and, thus, one array is needed to 
store fl. On the other hand, three arrays are necessary to store r\ at 
the past (n-1), present (n) and advance (n+1) stages. 

The actual vertical velocity itself, w, is not needed in the calcula- 
tion. However, when it is desired, it can be obtained from Equation (10), 
which is rewritten here in finite difference form. 

n+1 n-1 

"i.].^ H i.i a i.i.k-" I - 1 '"" ''i 2A ' t "i.[ 

M*i.k hi *M*vi.r x h i-i.rvM 

n+1 n-1 

+ H„ M + u * Vi.i ±hl 

2At i,j,k 2Ax 

+ v y , ^.i* " n J.T-i 

i.j.k 2Ay (36) 



The predictive equations for velocities are given by Equation (14) 
and Equation (15); in finite difference form, these equations are 



u n+1 u n+ \ - u 11 : 1 !!?: 1 . V. . . + V. , . .+ V. - . , . + V. . t „ 

l.l i.].k i,i i.i,k__. i.],k i-l,},k l-l, ]+l, k i.]+l,k 


2At " Iu i.j 4 




. n+1 n-1 
_„U r l,i " Vl.j . v \j.k-l " u i,i,k " a i,i,k 
" i.j Ax U. . (A a) 3 


+ U i.1.k + 1 „ 

" i.j.k 




(37) 


vf +1 v n+1 v - v n :K n -\ u. . . + u . . . . + 

i.l i,i. k i,] i.j.k __„, i.].k i,j-l.k 


u i+l,j-l,k + u i+l.j,k 


2At ^i.j 


4 


n. . - n. . , A v. . , , - vf 1 ! 1 . - v n ~\, 
cV i.l i .1-1 . v i.i.k-1 i,i. k 1,1, k 

b i.j Ay ■ V. . (As) 2 

1 »J 


'i.i.k+i 
-" '■ i. . . 

1.1 ,K 




(38) 



where, unless otherwise noted, all the superscripts are n and thus are 
dropped for convenience. The nonlinear inertia terms X. . . and Y. . . 
are 

A i,j,k " Ax 

, »W + 'l.lW.I+rt'.l+l ' <\.\+*l.H?l\ 

Ay 

n "l.i.k-l^i.i.k-l + "i-l.l.k-^ ' "l.i.k.l'ii.M * n i-i,i,k+l> 
+ U i.j 4Aa 

n+1 n+1 n-1 n-1 

+ [1 . 0pl)t . 1 °».i. fc -i-'H.i.^ M 'V..rYi -v.i 



(39a) 



Y. 



H i,i V i,i,k V i,i,k " H i,i-l V i.1-l,k V i,j-l,k 



i.j.k Ay 



(h i+i,i "Vi.i^ii.i.kVi.i.k- (h i,i +E i,i )u i^i.i.kVi.i,k 
+ v v i.i:k-i (n i,i.k-i * "u-i.k-i* - v i.i.k+i (i \i.k+i + n i,M,k+i } 

i.j 4Aa 

n+1 L n+1 n-1 n-1 

+ [l-(k-i)Ao] < j^ 4At 

(39b) 



Unless otherwise noted, all the superscripts are n and thus dropped, 
since the inertia terms are calculated at present time level n. In Equa- 
tions (39a, b) , n n and fl are presently obtained variables; since at each 
time level n and n are calculated first by Equations (34) and (35) respec- 
tively, u and v are followed by Equations (37) and (38). 

Whtin the calculations for velocities are to be perfprmed usincj-the 
predictive equations, (37) and (38), the variables u. * . and v. . . 
which jppeared respectively on the right-hand-side of 'Equations '(3?) and 
(38) should be rearranged, so that a form for velocities similar to Equa- 
tion (34) appears. In this form, all quantities on the right-hand-side are 
obtained based on specified h, calculated u and v at present and previous 
time levels, newly computed n and n at three time levels*' All the variables 
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and notations that appeared in Equations (37) through (39) have been 
defined and are shown in Figure 4. However, in Figure 4, the common 
superscript n and subscript k have been omitted for brevity. Figure 4 
also shows dearly the horizontal extent of variables involved in the cal- 
culation of n, u and v in the cross-hatched region, namely n. .# u, . . 
and v. . . (for k = 1 to 4), *'* '' J ' K 

Similar to the treatment of the momentum equations, the energy 
equation, Equation (16), can be written in finite difference form as 



„n+l_n+l „n-l_n-l 

i,3 i,],k i,] i,i, k _ 
2A t 

B h H ifj [ ^p »— + j— J, 1 

p T n _ T n+1 _ r 11 " 1 + T n 

. v x,j,k-l i.j,k i,i,k i,j,k+l _ n 
H. . (A a) 2 i,]\k (40) 

where the nonlinear term, R. . . , is given by 

R i.j,k = 4lx [u i+l,j (u ltl.j.k + Vl,J,k+l HT i+l.j.k + T i.j,kJ 



- U i.j (u i.i.k + u i,i.k + l )(T i.]\k + T i-l.j,k^ 



+ 4«iy tV i.j+l (v i.j+l,k + v i.j+l.k+l HT i,j+l,k + T i,].k ] 



- V i.j (v i.j.k +v i.J.k + l )(T i.J.k + T i.j-l.k )] 

H. . 
+ 2l7 2 i,j.k (T i.j,k +T i,j.k-1 ) " "i.j.k+l^i.j.k +T i,j,k+1 )] 

T T n+1 n-1 

+ [l-(k-*)Ac]. M'*' 1 " i.j.kH n i,j ' n l.) 



Likewise, unless indicated, all the superscripts are n and thus dropped 
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for clearity. Equation (40) , similar^tp Equations (37) and (38), can 
lead to a predictive equation for T* . ^. 

In the model, since the density is considered as a constant, the 
equation of state. Equations (17), is of no use, and the system is un- 
coupled. With that, the temperature field can be solved separately pro- 
vided that all results from the hydrodynamic model are available* That 
is, the spatial and temporal variation of temperature is solved after the 
spatial and temporal variation of velocity fields is known. n $Uernatiyely, 
Equation (40) can predict T after obtaining n , fl, u and v 
at each time cycle. The former involves two separated programs while 
the latter is a coupled program only. In coding the model above, both 
possibilities are taken into account by using flag statements (see Users 
Manual for Free-Surface Model, 1980). 

SOLUTION PROCEDURES 

Clearly, Equations (37) and (38) are for u and v at interior points; 
that is, they are not on the boundaries. Figure 4 implies why these 
formulae are not for the boundary points* Therefore, Figure 2 indicates 
that the normal velocity along the boundaries must be either specified or 
calculated by some other means- For solid lateral boundaries paral'el to 
x- or y-axis, the normal velocities are specified to be zero all the I'me. 
At the river mouth and discharge outlet, the normal components are given 
by the known flowrate and the average water aepth. At the seaward 
boundary, the imposing tide is specified along a parallel line at half-grid 
size away from the boundary. The difference of water elevations at 
half-grid points across the open boundary is used to calculate the normal 
velocity there. 

At the bottom, the no-slip condition requires that u = v = at all 
times; hence, there is no need for doing Equations (37) and (38) at K = 5. 
At the free surface or k = 1, the wind-produced shear stresses are pre- 
scribed as a function of space and time, and can be written as 

These conditions demand a modification of predictive formula for u n+ and 
v at k = 1. That is, the DuFort-Frankel format for second derivative, 
Equation (33), is replaced by (indices i and j are omitted and values of 
k are indicated) 

,3*u, n _ 2 n . n+1 . ,n-l . ^i.jV^ 
( 37 } k5sl ~ U7P [U 2 iU l * U 1 + ~ A^ 3 
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and a similar one for v-momentum equation* In this scheme, the free- 
surface .wind draa has been incorporated to drive the current through 
the u and v at the free surface. 

The boundary conditions on vertical velocity w are the kinetic con- 
dition 

• -U «t.-0 

and the rigid bottom condition 

w = at c = -1 

These conditions expressed in terms of modified vertical velocity a are 

ft = at a = or k = 1 

a = at a = -1 or k = 5 

As stated earlier, these conditions have been incorporated into the first 
integral form of the continuity equation; namely, the predictive equation 
for surface elevation. Equation (13) or Equation (14) with F = 0. Conse- 
quently, the finite difference formular for fi. Equation (35), is for k = 2, 
3 and 4, while Q k=1 and Q k=5 are set to zero at all times. 

The boundary conditions of energy equation are adiabatic on solid 
boundaries and known temperature or zero normal derivative of tempera- 
ture en open boundaries. It can be seen that if one neglects the lateral 
thermal diffusion terms, which are small in comparison with convection 
terms anyway, then no adiabatic condition is needed on the solid lateral 
wall; this is due to zero normal velocities on the solid lateral wall. How- 
ever, condition is still needed on the open boundaries- At the bottom, 
the heat flux is zero; while at the surface, the heat flux is proportional 
to (T - T ). These conditions are required when Equation (40) is used 
to calculate T at k = 1 and k = 4 respectively. 

As in most other hydrodynamic models for transient problems, the 
computation routine works step-by-step in time* This means that the 
computation proceeds through a sequence of time steps, each advancing 
the entire flow configuration through a small, but finite, increment of 
time, At- The results of the present and the last steps act as a basis 
for the calculation tc proceed to the next one, whereby the initial con- 
ditions can develop, within the limitations imposed by the boundary con- 
ditions, into the subsequent flow configurations. That is, provided that 
the values of dependent variables are known initially, the values at sub- 
sequent times are obtained by using the explicit scheme. The leap-frog 
finite difference formulae, Equations (34), (36), (37) and (40), predict 
surface elevation, n, and two horizontal velocities, u and v, and water 
temperature, T, at time level n+1. The vertical velocity w does not 
appear in these equations and thus, is left out until needed for picturing 
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the flow configuration* Each time cycle itself contains the following recur- 
sive processes. 

1. The countinuity equation in the form of Equation (34) is used for n 
through the mass conservation of the fluid £flumn at (i,j). At the 
same time, the modified vertical velocity Q is also calculated at 
k = 2, 3 and 4, since it is also based on the continuity equation but 
in a different form; namely, Equation (35)* 

2* The nonlinear inertia terms that appear in the horizontal momentum 
equations, (37) and (38), are considered as driving terms; that is, 
they are calculated from formulae (39a,b) by using known results 
at the present time (n) * This calculation usually takes a significant 
portion of the total computing time; thus, it is advised to drop this 
calculation whenever justified* For general tidal flows, the inertia 
effect is negligibly small in comparison with the Coriolis force. There- 
fore, if the Rossby Number, which is the ratio of inertial force to 
Coriolis force, is very close to zero, the calculation of X and Y terms 
can be skipped entirely. 

3. The horizontal momentum equations. Equations (37) and (38), are 
used to calculate u and v * Here, the specified boundary con- 
ditions, such as prescribed normal velocities, surface wind stresses 
and specified tide, come into play. 

4. The nonlinear terms, R, in the energy equation. Equation (40), is 
calculated by formula Equation (41). Unlike the nonlinear terms, 

X and Y of the momentum equations, R is to be included and calcu- 
lated under normal circumstances. 

5. Temperature at advance step, T # is calculated by Equation (40), 
and the adiabatic condition, known discharge temperature* given 
ambient temperature and surface heat transfer rate play a part in 
determining the temperature field. 

STABILITY 

The leap-frog (explicit) method has a limit on the size of the time 
step. Exceeding this limit makes the computation unstable. According to 
Platzman (I963), the maximum time step for an inviscid linear system, i.e. 
the system with viscous and nonlinear inertia effects neglected, the maxi- 
mum allowable time step is 

f - 2gH { 

max 2 (ax)* 

where H = maximum depth in the problem. If the Coriolis effect is 

also neglected, then 



it - Ax 

max 



J2^H 



max 

The leap-frog method with the DuFort-Frankel scheme for vertical 
diffusion terms has been adopted in the present model. As pointed out 
by Forsythe and Wasow (I960), this format is unconditionally stable for 
a pure diffusion model, a system in which the viscous diffusion is the 
only mechanism responsible for transport* Although, the present hydro- 
thermal model is considerably more complex than a pure diffusion model, 
the use of the DuFort-Frankel format is an important consideration for 
a shallow water system because the vertical diffusion criterion tends to 
become relatively more restrictive as the vertical dimension becomes 
smaller* 



28 



SECTION 5 
APPLICATIONS TO ANCLOTE ANCHORAGE 



INTRODUCTION 

The present model has been successfully applied in thermal disper- 
sion studies at the Anclote Anchorage site. The site, located on the 
west central coast of Florida near the town of Tarpon Springs, is 3 shal- 
low channel between the mainland and Anclote Key which separates the 
channel from the Gulf of Mexico, as shown in Figure 1. Anclote Anchorage 
is of interest since it is the receiving water body for the Anclote Power 
Plant cooling water discharge. The channel is relatively shallow with 
depths ranging from 0.3 to 3.6 m. Shallow regions of less than 0.6 meter 
comprise about 35% of the Anclote Anchorage area, which is approximately 
5 km in length and 6 km in average width. The principal driving mecha- 
nism for current circulation is tidal flux at the north and south entrances 
of the channel. The tide is predominantly semidiurnal with a mean range 
of two feet. Earlier measurements of temperature end salinity indicated 
the currents flow in and out through both entrances. However, the water 
exchange appears to be stronger in the south than in the north, or the 
currents generally flow in the north direction during flood tide and flow 
in the south direction during ebb tide. 

The Anclote Power Plant operated by the Florida Power Corporation 
has two 515 MW, oil-fired, electrical generating units. Cooling water is 
drawn from the Anclote River through a man-made channel. The six 
pumps delivering a total of 1,990,000 gpm (125.6 m 3 /sec) are designed to 
raise the water temperature of 2. 8°C above the ambient water temperature. 
The heated water is discharged back into the Anclote b?< : n through the 
discharge channel with dredged submarine extension. T v* designed 
total flowrate is approximately S3 times th*^ long-tenn average flowrate 
of the Anclote River. At present, only Unit 1 is operative, while Unit 
2 is awaiting permission. The present flowrate, therefore, is 995,000 
gpm (62.7 m 3 /sec). 

The model as applied to the Anclote Anchorage shows its capacity 
for considering the effects of geometry and bathymetry, spatio-temporal 
variation of the free surface, various boundary conditions, including tides 
of different phase and range, surface heat transfer based on equilibrium 
temperature concept, and changing meteorological conditions. In addition, 
turbulence has been considered by using the eddy transport approxima- 
tions* 
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CHOICE OF DOMAIN AND GRID SYSTEM 

In this study, the Anclote Anchorage is designed to contain the area 
of water bounded by the Anclote Keys and sand barriers on the west, 
mainland on the east, and two imaginary E-W lines drawn from the northern 
tip of the barriers and the south end of Anclote Key to the mainland. 
These two imaginary boundaries are considered to be so far that the 
thermal plume will not reach them even at low tidal stage. Its location 
and schematization are shown in Figure 7. This water Is open to tides 
from the Gulf of Mexico at both ends. The Anclote River, where the 
intake of cooling water is located, is also included since the recirculation 
of discharc 'low is of concern to the power company. 

The cdopt^o horizontal grid layout and its index are also shown in 
Figure 7. The grid work is allowed to orient away from north-south, 
east-west system, but in general, the y-axis of the grid system aligns 
with south-north. Thus, the subscript i increases eastwardly while j 
increases northwardly. The z-axis is chosen upward from mean water 
surface while the subscript k increases downward from the water surface. 

The selection of the grid size is governed by several constraints, 
if the grid size is too large, the approximation of the channel in the 
system will be inaccurate, and at certain size, the computation will become 
meaningless. However, decreasing the grid size will iead to a considerable 
increase in the computer time since the computations must then be made 
on more points. In addition , the time step will decrease because the 
dispersive properties of the computational method are related to the ratio 
of the time step to the spatial grid size. After some numerical experiments, 
the model on a 16 x 14 x 5 grid with grid size A x = A y = 416. 75 m and 
A a - 0. 25 is considered a good compromise between the resolution desired 
for the region near the discharge and the limitation of the computer 
(UNIVAC 1108) at the University of Miami. Also, care has been taken 
to hzve the intake and outlet at grid points. The velocities at these 
grid points are specified such that the flowrate and direction can be 
easily represented. 

SUMMARY OF DATA 

Field measurements of current velocities and water temperature have 
been made for model calibration as well as verification purposes for the 
hydrothermal prediction. Two field trips were carried out in the summer 
and winter respectively. The procedures and data results will be discussed 
briefly below. 

June 1978 Data Acquisition 

On June 19 and 20, a team was sent out to the field for measurement 
of current velocity and temperature. In coordination, the infrared (1R) 
scanner data was obtained by flights over the channel. The current 
measurement and temperature readings were obtained at 9 points; 4 
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points were located at the north end, four others at the south end, and 
one in the middle of the channel, as shown in Figure 8. At each point, 
the measurements were done at different depths, namely the water surface 
and some depths below the water surface. 

The in-situ measurements were carried out by personnel on three 
boats. The Barnes PRT-5 Radiometer was used to measure the surface 
temperature, and a thermister (#9) with floating mechanism was used to 
measure the temperature of the water surface. Temperature profiles are 
measured by lowering another thennister (#4) into the water in increments 
of 2 to i feet. Current velocities are measured with a Bendix Model No. 
665 w /readout current meter which reads the magnitude and direction of 
the current. Figures 9 through 12 show the current velocities at 4 tidal 
stages, which give a picture of the current patterns in the channel, 
expecially at the boundaries. One may notice that the current is stronger 
at the south boundary than at the north boundary. Both boundaries can 
exchange water with the Gulf of Mexico. 

The flights were made at 609.6 meter or 2000 feet altitude. The 
black body of IR scanning window was set at different ranges for each 
flight set. For the first flight, the black body range was 74° to 98°F, 
the second was 75° to 99°F, the third was 73° to 102°F and the fourth 
was 80° to 104°F. All have a full black body range of 24°F; it provides 
a satisfactory resolution of 4°F for each of the six color? between white 
and black. A finer resolution was also obtained which reduced the color 
band to 0. 74°C (1.33°F) temperature spread. At 2000 feet altitude, the 
scanning width is 941.8 meters or 3090 feet. With this scanning width, 
the whole stretch of the channel was covered by ten east-west flights 
as labeled in Figure 13. 

The in-situ measurement of water surface temperature at the time 
when the airborne IR data was undertaken provides a calibration of IR 
temperature. Figures 14-17 show the measured temperature at their 
points at four tidal stages of this mission. Since the channel is quite 
shallow and only the surface temperature is of primary concern in this 
study, only the surface temperature is presented. 

January 1979 Data Acquisition 

Field measurements by boat and by IR scanning were schedualed 
on January 30 and 31. On previous experience, the current measurement 
and temperature reading were carried at different points from the pre- 
vious mission. Figure 18 shows the location of all points where the data 
were collected by three boats. The route for boat #1 was along points 
4, J, 2, ^ and S, and for boat #2, the route was along points 8, 7, 6, 
10, and 9. Meanwhile, boat #3 was working on the region near the coast 
between outlet and intake. The water temperature and current velocity 
were measured by the same instruments used in the previous data acqui- 
sition. The measurements were taken at the surface and successive 
depths of 3 ft (1 m) intervals. The flights were coordinated at the same 
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time as for IR scanning, while the boats were collecting the data. These 
ten flight lines, shown in Figure 13, were made at 609*6 meters or 2000 
feet altitude. The black body of IR scanning range was set at 44-80°F 
so each of the six color bands between white and black would represent 
a 6°F interval. 

On January 30, the weather permitted both morning and afternoon 
data collection. However, on January 31, the conditions was so bad that 
the mission had to be postponed to February 1 afternoon. Figures 19 
through 21 show the current velocities at three tidal stages; namely, 
flood tide, ebb tide and high tide respectively. Figures 22 through 24 
show the surface temperature distribution at corresponding stages. To 
further show the measured temperature fields, isotherms were interpolated 
from these in-situ data and presented in Figure 25 through 27 correspond- 
ing to each tidal stage respectively. 

In these two data acquisitions, although a plan of synchronized 
measurements of current at tide changes was carried out in order to 
provice current data at slack, flood and ebb, the effort was not so 
successful due to technical difficulties in ootaining reading and other un- 
expected circumstances. Therefore, the current data could only be of 
use as reference. 

CALCULATION OF INPUTS 

The important input parameters and some specification of boundary 
conditions, such as intake and discharge velocities, discharge temperature, 
tidal condition, river flow, surface heat transfer, and wind stress, will 
be presented in this section. 

1. Time step, DT 

In order to determine the time step, DT, the stability criterion 
has to be followed: 

Dx 41760 

DT < ;/2gH™ = j2x980x360 = 50 sec 

M * max " 

About 1/3 of this value is reasonably safe to use. 
Here, we use DT =15 sec. 

2. Vertical eddy viscosity, A 

The vertical eddy viscosity is estimated by means of the formula 

A v = CH 4/3 

where H is the local depth and C is an empirical constant. This 
type of formula for horizontal turbulent diffusion was originally 
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suggested by Richardson (1926). Since that time, the 4/3 power has 
been substantiated by a considerable amount of empirical evidence and 
theoretical analysis. The so-called constant C may vary with the in- 
tensity of the turbulence and is not well established. In this study, 
H is used as maximum depth and C is 0.002, so that we have 

A y = .002x(360) 4/3 = 6 cm 2 /sec 

For shallow, well-mixed tidal water, about three times the calculated 
value was found suitable. Here, we use A =20 cm 2 /sec. 

3. Horizontal eddy diffusivity, B. 

The horizontal eddy diffusivity is calculated by the same formula 
as mentioned above. However, the maximum depth is replaced by the 
maximum length of the domain, which is 6 km in this study. So we 
have 

B h = .002x(600,000) 4/3 = 100,000 cm 2 /sec 

4. Vertical eddy diffusivity, B 

In this study, the turbulent Prandtl Number is assumed as 1. 
Thus, the vertical eddy diffusivity is equal to the vertical eddy vis- 
cosity or B = A =20 cm 2 /sec. 

* V v 

5. Surface heat transfer coefficient, K 

s 

The procedures for K calculation are as follows: 

a - T . = T - (14.55 + C.114T )(l-f) - [(2.5 + 0. 007T )(1-f)] 3 
g 3 a a 

where T . = dewpoint temperature, °F 
T - air temperature, °F 
f = relative humidity in fraction of unit 

b. 8 = 0.255 - 0.0085T + 0.0002G4T 2 

ave ave 

where T ave = (T s + T rf ) /2, and S is an intermediate step 
T = ambient surface temperature, °F 

c. f(u) s 70 + 0. 7u 2 

where u = wind speed, tnph 

d. K s = 15.7 + (0 + 0. 26)f(u) 

where K $ is the surface heat transfer coefficient in BTU/ft 2 day. 
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The values of T , f, u, T needed for this calculation are read 
from the climato&gical data s . 

6. Wind drag coefficient, C rf 

The wind stresses on the free surface are introduced into the 
model as 

A A 

v 3ui _ _y 3_V| 

pTi 3a'o=0 " \ pH 3a'o=0 ~ T y 

The subscripts x and y indicate the shear stress acting in the x and 
y direction respectively. The relation of these stresses to the wind 
speed at a certain height is very difficult to determine theoretically 
and its value is usually based on semi-empirical formulae. The well 
known form of the relationship between shear stress x and wind 
speed U (usually measured at a height of ten meters) is 

where p is the air density and C. is a dimensionless drag coefficient. 
in this study, the drag coefficient formulae obtained by Wu (1969) is 
adopted, and its value was given by 

C , = 0.00125 U 2 for U < 1 m/sec 
d 

= 0.0005/0" for 1 < U < 15 m/sec 
= 0.0026 for U > 15 m/sec 

7. Intake and discharge velocities 

The intake and discharge velocities are calculated according to 
the discharge flowrate from power plant data, the grid size and the 
average depth at the intake and discharge outlet. The procedures 
are shown as follows: 

a. Flowrate = 955,000 gpm (from power plant physical data) 

= 62.8 m 3 /sec 

= 62.8 x 10 b cm 3 /sec 

b. Both intake and discharge channels are at 45° from N, 
therefore, 

31.4 x 10 cm 3 /sec is crossing the Ax and Ay at the point of 
intake and discharge. 

c. The average depth at intake and discharge outlet is approxi- 
mately 4 1 or 122 cm, and the width is ix = Ay = 41760 cm; 
so the cross-sectional area is 41760 x 122 cm 2 . 
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d. The average velocities: 



U = V 
ave ave 



31.4 x 10 v 



41760 x 122 
The velocity profiles are assumed as shown. 



= 6.163 cm /sec 



7.0 cm /sec 




T' 
i J 

k 3 

4 

5 

f. To allow for channel storage during tide change we assume the 
intake and discharge velocities to be sinusoidal, i.e. 

Intake: V3(14, 4, k) = 7 - 3 x cos[j~ (EST - 7.625)] 

U3(15, 3, k) sV3{14, 4, k) for k = 1, 2, 3, 4 

Discharge: V3(14, 8, k) = 7 - 3 x cosl-—^ (EST - 7.5)] 

U3(15, 8, k) = -V3(14, 8, k) for k = 1, 2, 3, 4 

where 7.625 and 7.5 are taken to be the values of the phase shift 
which takes into account the time to travel from the south end of 
Anclote Key to the concerned point. 

8. Tidal condition on June 19, 1978 

Simulated diurnal tide is shown in Figure 28, where 

a. Period = 12.5 hr 

b. Stage = short term average sea level - MSL = 48 cm 

c. Amplitude = i short term average tide range = 65 cm 

d. Time shift = 7*125 hr 

i.e. at 7,125 am, June 19, 1978, the tide at the south end of 
Anclote Key was zero. 

e. W - E lapse = 0.014 hr/DX 



Wave propagation speed C =J 2gh -J 2x980x360 = 850 cm /sec 
(H = 360 cm is the maximum depth of the Anchorage.) 
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The time needed to travel one grid distance is 
DX 41760 



350 



= 50 sec = 0.014 hr 



We use 0.014 hr per DX for phase shift in W - E direction and 
the imposing tide at the south entrance is 

n s = 48 + 65 sin[y|^ (EST - 7.125 - 0.014(1 - 1)] 

I = grid no. in W - E direction. 

f. S - N lapse = 0.15 hr 

Distance from souift entrance to north entrance is about 543000 cm. 

Time for wave to travel this distance is 5a30QO = 0.18 hr. 

We take 0.15 hr as phase difference betwee?! the south and the 

north boundaries; there, the imposing tide at the north entrance 

is 

n =48 + 65 sin[~c ( EST - 7.125 - 0.15 - 0.014(1 - 1)] 

9. Tidal condition on January 30, 1979 

Simulated diurnal tide is shown in Figure 29. The calculation 
procedures are of the same as summer tidal conditions. 

a. Period = 12.0 hr 

b. Stage = 36. 6 cm 

c. Amplitude = i short-term average tide range = 42.7 cm 

d. Time shift = 10 hr 

i.e. at 10 am, January 30, 1979, the tide at the south end of 
Anclote Key was zero. 

e. W - E lapse = 0.015 hr/DX 

This value is slightly higher than the summer case since the 
maximum water depth in the winter is less than the maximum water 
depth in the summer. The imposing tide at the south entrance is 

n s = 36.6 + 42.7 sin[y| (EST - 10 - 0.015(1 - 1)] 

f. S - N lapse =0.2 hr/DX 

This is the time for wave to travel from south entrance to north 
entrance. So the imposing tide at the north entrance is 

n = 36,6 + 42.7 sin[~? {EST - 10 - 0.2 - 0.015(1 - 1)] 
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10. Anclote River flowrate and temperature 

a. The distance traveled from South Anclote Key to Tarpon Springs 
is 20 DX. We estimate a time lapse of 0.5 hr to account for the 
retardation due to buffering effect of river storage and Anclote 
River's natural outflow. 

b* The average current is estimated to be 20 cm /sec, 
therefore, we take 

U3(16, 1, k) = 20 COSty^j (EST - 7.625)] 

V3(15, 1, k) = -20 cost—^ (EST - 7.625)] 

for k = 1, 2, 3, 4 

c* The surface elevation at Tarpon Springs is to be calculated. 

d. To be in accordance with given velocities at Tarpon Springs, 
the temperature there is also assigned and its value has a 
24 hr period instead of 12.5 hr. This temperature on June 
19, I978 is 

T3(15, 1, k) = 26.9 + 0.5 sin[|j (EST - 12)] 

while on January 30, 1 979 the temperature was assumed as 

T 3(15, 1, k) = 11.9 + 0.5 sin[|| (EST - 12)] 

where the 12 hr shift is to make the peak temperature occur 
at 1800. Thus, the water in and out at Tarpon Springs has a 
temperature ranging from 26.4 (before dawn) to 27.4 (late after- 
noon) in the summer and ranging from 11.4 (before dawn) to 
12.4 (late afternoon) in the winter. 

11. Discharge temperature 

On June 1^20, 1978, the recorded discharge temperature at 
daytime was in the range of 29.3 - 30. 3°C, while on January 30- 
February 1, 1979, the temperature range was 16. 4 - 15.2°C. To 
account for the further drop of discharge temperature due to cooler 
ambient temperature at nighttime, we assumed a sinusoidal variation 
of discharge temperature with diurnal period. 

a. Discharge temperature is estimated for June 1978 

T3(14, 8, k) = 29.4 + 0.4 sin(|| (EST - 12)] 

Therefore, the highest discharge temperature of 30.3°C happens 
at 6 pm and the lowest 29. 3°C at 6 am. 
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b. Discharge temperature is estimated for January 1979 

T3(H, 8, k) = 15.8 + 0.6 sin(~ (EST - 12}] 

The highest discharge temperature is 16. 4°C at 6 pm and the 
lowest discharge temperature is 15. 2°C at 6 am. 

12. The Gulf temperature 

The Culi water outside the Anclote Anchorage as well as the 
atmosphere is sink to the heat disposal from the power plant; 
therefore, the boundary conditions on temperature at the north and 
south entrance are not considered as adiabatic as in normal case of 
far-field thermal pollution problem. Instead, we specify the outside- 
Anchorsge ambient temperatures. Again, they are 24 hr periodic 
and their values should be in accordance with the measured tempera- 
ture in the same neighborhood. Here in compliance with measured 
data, we use 

T ab = 27 '° + °* 2 sinf l5 (EST " 12)I 
for the summer simulation during June 19-20 I978 and use 

T ab = Um8 + °* 4 sin[ 4? (EST ~ U)] 
for the winter simulation during January 30-February 1, I979. 

RESULTS 

Due to the fact that the archival data for Anclote site is inadequate, 
the following inputs are used for the computer runs: data from NOAA 
tide table, solar radiation, wind and the power plant operating conditions. 
Tables 1 and 2 show the operating conditions of Anclote Power Plant 
during June 19-20, 1978 and January 30-February 1, I979 respectively* 
The isotherms obtained from the IR data from each flight are interpolated 
by hand from mosaic digicolor films and then plotted by the computer. 
These isotherm plots are presented in this section for easy comparison 
with the predicted isotherms. The average deviation of the calculated 
temperatures from IR temperatures is also indicated. This deviation is 
calculated by simply averaging the temperature differences between the 
measured and calculated temperatures at each point of the domain. The 
average deviation is given by 

6 2 = S[TB(i, j, 1) - TIRO, j)I 2 /I (i, j) 

where TB is the calculated temperature, TIR is the IR temperature and 
£(i, j) is the number of surface half-grid points in the detain. 
U 

The TIR and TB isotherms are then compared to assess the accuracy 
of the modei in predicting the dispersion of the waste heat. It is recog- 
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nized that the principal factors affecting the flow pattern and the shape 
of the thermal plume are the tide driving through the north and south 
boundaries, wind effects on the water surface, bottom topography, heat 
transfer through the air-water interface, and the intensity of eddy vis- 
cosity and diffusivity* The effects of each of these factors on the flow 
pattern have been discussed in detail by Lee et al. (1978a, b) . These 
effects are important in understanding the numerical behavior of the model, 
and play a leading role in proving the capabilities of the model. In this 
report, however, only the results of the verification runs are presented. 
The results of the summer and winter verifications are discussed in the 
next two sections. 

Summer Results 

The figures from the summer simulation show the hydrothermal dis- 
persion of waste heat under the conditions of June 19-20, 1978, as de- 
scribed in the previous section. The simulation run started at 0400 EST, 
June 18, with initial conditions of zero velocity field, equilibrium water 
level and constant ambient temperature (so-called "cool start 11 ). The tide 
and heated discharge are then imposed. After 20 hours of "warm up," 
the thermal plume can be seen to develop. Experience indicates that the 
initial condition is not important, as its effects die out in the first one 
or two hours. However, it is convenient to begin a simulation through a 
"cool start." 

Figure 30 shows the surface flow pattern at 1030 EST, June 20, 
which corresponds to the high tide at south Anclote Key. The current 
velocity is relatively small and the mainstream flows in the north direc- 
tion. A recirculation can be seen to occur as a result of the flow of 
water from the Anchorage into the Anclote River. The net effect is the 
flow of part of the discharge into the river entrance where the intake 
structure is located. 

Figure 31 shows the resultant velocity at high tide. The result of 
the velocity components, u and w, is plotted on the vertical cross sections 
at J = 4, 8 and 12. It is believed that the currents at these sections are 
most affected by the plan- form configuration, bathymetry and other fea- 
tures such as the river mouth, discharge and tidal boundary etc. The 
u-w velocity profiles show clearly the effects of bottom friction in retard- 
ing the flow. Figure 32 shows the v-w resultant velocity on the vertical 
cross sections at 1 = 4, 8 and 12- In these two figures, the vertical 
velocity component w has been exaggerated to make the vertical circula- 
tion of the current detectibie. One may notice that, since the velocity 
vector at the top level (c = 0) is plotted right on the water surface, 
joining the tails of these velocity vectors would show the free-surface 
profile. 

Figure 33 shows the isotherms plotted from the IR data. The com- 
puter-plotted points were obtained through visual interpolation from the 
mosaic IR image, as discussed earlier. Due to the difficulties of visual 
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interpolations and large grid size, the computer-plotted isotherms un- 
avoidably differ somehow from the real isotherms in the digicolor film. 
This figure shows that the heated water is recirculated back to the intake 
through the Anclote River. This will cause an increase in the input 
temperature of the condenser cooling water, thus reducing the efficiency 
of the power plant. 

Figure 32 shows the computed isotherms at the same tidal stage* 
Comparing this figure with Figure 33, the good agreement of these 
temperature fields can be easily seen. The average deviation of the 
simulated temperature from IR temperature is only 0. 359°C. 

Figures 30 through 34 form the first set of results of the 
summer simulation. This set of plotted results, which includes the 
surface velocity, u-w velocity, v-w velocity, IR temperature and the 
calculated temperature, coincides with the high tide at south Anclote 
Keys, Figures 35 through 39 is the second set of results of the same 
simulation run, at 1430 EST, June 20, or at a time corresponding to 
maximum ebb tide at south Anclote Keys. Therefore, this set shows 
the current and temperature fields of the Anchorage four hours after 
the first set of results was taken. It is also seen that the Anclote 
River empties into the Anclote basin while the basin drains northward 
and southward. There are strong ebb flows occurring at both the north 
and south end of the Anclote Keys barriers. The dividing line lies along 
the j = 10 grid line. On the northern side of this line, the water flows 
north, while on the southern side, it flows south. At this ridge line, 
the water has minimal transverse motion. The legation of this line varies 
with time as the tides from both ends are not in phase. In reality, most 
of the flow in the Anclote Anchorage comes primarily from the south end; 
therefore, the ridge line is observed to exist at a region close to the 
northern end of the Anchorage. Also, this line is observed to move 
southward during ebb tide and northward during flood tide. The flow 
pattern corresponds with the observed flowfield at a similar tidal stage. 

Figure 36 shows a strong u-velocity component at the river mouth, 
while Figure 37 shows a strong v-velocity component near the south 
boundary. Figures 38 and 39 are two corresponding isotherm plots at 
this tidal stage. It can be seen that the thermal plume moves according 
to the stage of the tidal cycles, as would be expected. The 27. 50°C 
isotherm from the simulated results covers a larger area than that of the 
IR plot. Generally speaking, the simulated results are in yood agreement 
with the IR data. 

Figures 40 through 45 is the third set of results of the simulation 
run, but at 1730 EST, June 20, which corresponds to low tide at the 
south Anclote Keys. The current field is seen to be generally quite small 
except at Anclote River region. The river is still undergoinq its outflow- 
ing process and continues to empty its tidal storage into the basin, while 
the open sea flows into the Anchorage through both open boundaries. 
It has to be mentioned here that the lower part of the 27. 50°C isotherm 
from the IR data is not very visible in the IR digicolor film. As a result, 
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the contour shown in Figure 43 is estimated by extrapolation from the 
IR data. The calculated isotherm at this temperature, as shown in 
Figure 44, has a peculiar shape at the same region* This may be due 
to the fact that this region, north of the dredged ship channel leading 
to Tarpon Springs, is relatively shallow and becomes especially so during 
low tidal stage. The sum effect is that the heat transport, due to con- 
vection, becomes smalK This shallow region can be easily seen from 
Figure 42 along the I = 8 cross section at a location about J = 4, 5 and 
6. If this isotherm is ignored, then the remaining isotherms, shown in 
Figures 43 and 44, have the same sort of tongue-shape profiles, show- 
ing some degrees of good agreement. 

Figures 45 through 49 is the last set of results of the summer 
simulation. The time is 2030 EST at the same day and corresponds to 
flood tide at the south end of the Anchorage. The incoming tide from the 
south end drives northward into Anclote basin. The dividing line dis- 
cussed earlier, does not appear, as the northbound tidal current changes 
course into the northeast direction. Both Figures 48 and 49 show that 
the thermal plume has been pushed back, and the isotherms become more 
compact* These features could be explained by considering that the 
changing tide retards the convective transport of the thermal plume. 
The 27,50 and 28, 25°C isotherms in Figure 49 are seen to be pushed 
toward the northeast direction rather than being pushed toward the east 
coast, as shown by the corresponding IR isotherms, Figure 48. This 
could be because the north and south tidal conditions are no longer con- 
sistent with the real tidal conditions after such a long time of simulation 
(64.5 hours). At this stage, it seems that the northern tide could in- 
duce a stronger current to push the plume toward the southeast direction. 
If correct tidal data had been provided, the results might have improved. 

Winter Results 

The winter simulation run started at 2200 EST, January 29, 1979. 
The meteorological data and input conditions are as discussed in the pre- 
vious section. Here, we present only two sets of results corresponding 
to the successive flood tidal and maximum ebb tidal stages at the Anchorage 
on January 30, 1979, These tidal stages were recorded during the first 
and second data acquisition missions which took place en January 30, 1979, 
at about 1100 EST and 1600 EST respectively. It should be pointed out 
that, due to sudden change in weather condition, the second and third 
missions were about 48 hours apart. In between these missions, the sea 
had become very rough under very stormy conditions. The effect of the 
stormy weather can be clearly seen by comparing the two IR temperature 
fields, Figures 58 and 60; the former was taken at about 1600 EST, 
January 30, while the latter was taken at about 1600 EST, Feburary 1. 
There is therefore a 48 hour time lag. Figure 58 shows the thermal 
plume at ebb tide of a typical winter day. This is a great contrast to 
Figure 60 which shows the plume during high tide but after the in- 
fluence of the stormy weather. It is to be noted that the isotherms 
shown in Figure 60 are not of the same values as those shown in Figure 
58, In fact, there was a 4°C drop in temperature, and a 2. 3°C drop in 
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the discharged water temperature. Particularly, the wind condition had 
gone through a very substantial change during this 48 hour time gap. 
This wind speed, shown in Figure 60, was still 50% stronger than that 
shown in Figure 58. Therefore, it was no surprise to find that the 
thermal plume, shown in Figure 60, was so compact and shrunken. 
The windy condition and the relatively low ambient temperature were 
the apparent causes of this much dwarfed plume. Observing the com- 
pactness of the neighboring isotherms and judging the resolution power 
of the present grid system, one is tempted to suspect the ability of 
modeling under this kind of extremely severe meteorological condition. 
Experience shows that it is difficult to obtain such a compact thermal 
plume through simulations of such adverse natural conditions. Consider- 
ing the roughness of the grid system, this is actually not a surprise. 
Therefore, for the winter simulation only two sets of results, namely, 
1100 EST and 1600 EST on January 30, are presented. 

Figures 50 through 54 is the first set of results of 1100 EST, 
January 30, 1979, which corresponds to the flood tidal stage at south 
Anclote Key. Figure 50 shows the surface-flow pattern from calculated 
results. The incoming tide from the south entrance drives directly into 
the Anclote River. There is a dividing line across the Anchorage at 
this tidal stage. This ridge line ic located along the J =9 grid !ine, 
at which the incoming tides from botr north and south meet. This can 
also be seen from the in-situ measurement, as shown in Figure 19* The 
discharge at the outlet is pushed head-on by this flew. As a result, 
the thermal plume is squeezed by this flow. The isotherm plots in 
Figures 53 and 54 show this effect very clearly. Although Figures 53 
and 54 show similar tendencies, the calculated isotherms seem not to be 
in good agreement with the IR isotherms; particularly so for the 15. 40°C 
isotherms. However, the calculated plume areas are close to those ob- 
tained from the IR data. Both temperature fields show the dispersion 
of the thermal plume along the shore. The recirculation of heated water 
is clearly indicated in both figures, showing that a recirculation actually 
occurs at this tidal stage under the meteorological influence. 

Figures 55 through 59 show the second set of results of the winter 
simulation runs. This set of plotted results shows the current and 
temperature fields at 1600 EST, January 30, 1979, during ebb tidal con- 
ditions at the Anchorage. Figure 55 shows a strong current driving from 
northeast to southwest along the channel. No ridge line is observable, 
and the discharged water is not entrained by this main current; in fact, 
most of the discharged water flows directly through the river mouth into 
the intake of the power plant. Thus, serious recirculation happens at 
this time; this is also shown in Figure 58 of the IR isotherms and in 
Figure 59 of the calculated isotherms. The calculated current pattern 
shown in Figure 55 is in good agreement with that of the in-situ mea- 
surements shown in Figure 20. A similar agreement can be seen between 
the calculated thermal plume shown in Figure 59 with that of the IR data 
shown in Figure 58. It can also be seen that all the calculated isotherms 
cover a slightly larger area than the corresponding IR isotherms; this 
difference is shown to be insignificant by a low deviation of 0.65°C above 
the IR measured temperature. 
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TABLE1 CL I MATIC DATA FOR SUMMER RUN AT ANCLOTE ANCHORAGE 



Date 


EST 


(1) 
T air(°C) 


(2) 

Humidity 


(3) 
Wind 
Soeed 


(4) 
Wind 
Direction 


(5) 
Solar 
Radiation 


(6) 

T 
surf 


6 /IS 


04 


23.3 


.90 


357.6 


50. 


0.0 


26.4 




05 


23.3 


.90 


268.2 


50. 


0.0 


26.4 




06 


23.3 


.90 


312.9 


50. 


0.0 


26.5 




07 


23.3 


.90 


312.9 


50. 


0.05 


26.5 




08 


25.0 


.84 


312.9 


50. 


0.40 


26.7 




09 


26.7 


.77 


447.0 


70. 


0.75 


26.8 




10 


28.3 


.70 


536.4 


80. 


1.05 


26.9 




11 


30. 


.63 


536.4 


80. 


1.40 


27.0 




12 


31.1 


.60 


536.4 


70. 


1.6*0 


27.0 




13 


32.2 


.56 


536.4 


90. 


1.70 


27, 1 




14 


31.6 


.57 


491.7 


80. 


1.60 


27.1 




15 


31.1 


. 58 


536.4 


80. 


1.50 


27.2 




16 


30.6 


.59 


581.1 


90. 


1.30 


27.2 




17 


30.0 


.61 


581.1 


80. 


1.10 


27.2 




18 


29.4 


.64 


625.8 


80. 


o;7o 


27, 1 




19 


28.3 


.67 


581.1 


SO. 


0.30 


27.0 




20 


28.3 


.70 


447,0 


SO. 


0.05 


26.9 




21 


27.2 


.73 


402.3 


80. 


0.0 


26.9 




22 


26.1 


.77 


312.9 


SO. 


0.0 


26.3 




23 


25.0 


.75 


402.3 


30. 


0.0 


26.3 




24 


23.9 


. 74 


312.9 


30. 


0.0 


26.8 


6/19 


01 


22.2 


.73 


223.5 


SO. 


0.0 


26.8 




02 


21.7 


.76 


223.5 


60. 


0.0 


26.8 




03 


21.1 


.80 


312.9 


40. 


0.0 


26.8 



16 









TABLE 1 


(Continued) 






Date 


EST 


(1) 
T air(°C) 


(2) 

Humidity 


(3) 
Wind 
Soeed 


(4) 
Wind 
Direction 


(3? 

Solar 
Radiation 


(6) 

T 
sun 


6/19 


04 


21.1 


.34 


ol2.9 


50. 


0.0 


26.8 




05 


21.7 


.82 


312.9 


50. 


0.0 


26.9 




06 


22.2 


.30 


312.9 


70. 


0.0 


26.9 




07 


22.3 


.79 


312.9 


70. 


0.10 


27.0 




08 


23.9 


.73 


357.6 


70. 


0.25 


27.0 




09 


25.0 


.67 


402.3 


SO. 


0.40 


27.0 




10 


26.7 


.60 


447.0 


80. 


0.60 


27.0 




11 


27.3 


. 56 


491.7 


90. 


1.25 


27.0 




12 


29.4 


.53 


536.4 


SO. . 


1.25 


27.0 




13 


30.6 


.50 


447.0 


80. 


0.60 


27.0 




14 


30. 


.53 


312.9 


90. 


0.30 


26.9 




15 


29.4 


.57 


268.2 


SO. 


1.20 


26.9 




16 


2S.3 


.61 


447.0 


350. 


0.70. 


26.9 




17 


28.3 


.65 


536. 4 


350. 


0. 50 


26.9 




18 


27.2 


.69 


581.1 


360. 


0.15 


26.9 




19 


26.1 


.74 


581.1 


350. 


0.15 


26.8 




20 


25.6 


.76 


357.6 


60. 


0.05 


26.8 




21 


25.0 


.78 


402.3 


120. 


0.0 


26.7 




22 


25.0 


.79 


357. 6 


110. 


0.0 


26.7 




23 


25.0 


.79 


312.9 


110. 


0.0 


26.6 




24 


24.4 


.79 


357.6 


110. 


0.0 


26.5 
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TABLE 


1 (Contin 


uedl 






Date 


EST 


(1) 
T air(°C) 


(2) 
Humidity 


(3) 
Wind 
Speed 


(4) 
Wind 
Direction 


(5) 
Solar 
Radiation 


(6) 

T 
surf 


0/20 


01 


24.4 


.79 


312.9 


100. 


0.0 


26.4 




02 


23.9 


• 30 


268.2 


100. 


0.0 


26.5 




03 


24.4 


.31 


357.6 


110. 


0.0 


26.6 




04 


25.0 


.32 


357.6 


110. 


0.0 


26,6 




05 


25.0 


,S2 


402.3 


110. 


0.0 


26.7 




06 


25.0 


.32 


357.6 


no. 


o.a 


26.7 




07 


25.0 


.82 


357.. 6 


90. 


0.05 


26.7 




OS 


26. 1 


.73 


357.6 


' 90. 


0.20 


26.3 




09 


27.2 


.74 


491.7 


110. 


0.60 


26.8 




10 


2S.3 


.70 


581.1 


110. 


0.30 


27.0 




11 


23.9 


.66 


447.0 


110. 


0.40 


27.0 




12 


23. 4 


.62 


268.2 


110. 


0.40 


27.0 




13 


30.0 


.59 


357.6 


90. 


0.60 * 


27.0 




14 


30. 6 


.53 


402.3 


90. 


0. 50 


27. 1 




15 


30.6 


.57 


402.3 


110. 


0. 55 


27.0 




16 


30.0 


.57 


268.2 


100. 


0.40 


27. 1 




17 


30.0 


.56 


402.3 


100, 


0.30 


27 2 




IS 


29.4 


.55 


312.9 


110. 


0.20 


27. 1 




19 


23.8 


.53 


22o. 5 


90. 


0. 15 


27.0 




20 


23,0 


. 50 


223. 5 


90. 


0.05 


27.0 
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T 


ABLE 2 


CLIMATIC 


DATA FOR WINTER RUN AT ANCLOTE ANCHORAGE 






Date 


EST 


(1) 
T air(°c) 


(2) 

Humidity 


(3) 
Wind 
Speed 


(4) 
Wind 

Direction 


(5) 
Solar 
Radiation 


(6) 

T 
surf 




1/29 


22 


a.s 


♦ 98 


350 


320 


0.0 


11.8 






23 


8.8 


.93 


350 


320 


0.0 


11.3 






24 


8.3 


.89 


350 


320 


0.0 


11.8 




1/30 


1 


8.3 


.36 


350 


330 


0.0 


11.7 






2 


7.7 


.87 


350 


340 


0.0 


11.7 






3 


7.2 


.88 


325 


350 


0.0 


11.7 






4 


6.7 


.89 


300 


360 


0.0 


11.7 






5 


6.1 


.39 


230 


10 


0.0 


11.7 






6 


5.5 


.89 


270 


20 


0.0 


11.7 






7 


5.0 


,S9 


250 


40 


0.0 


11.7 






8 


6.7 


.75 


230 


60 


0.05 


11.7 






9 


3.3 


.63 


310 


80 


0.25 


11.3 






10 


10.0 


.54 


350 


100 


0.45 


11.9 






11 


11.6 


.50 


375 


140 


0.75 


12,0 






12 


13.3 


..47 


400 


180 


0.95 


12.1 






13 


15.5 


.43 


400 


220 


1.15 


12.2 






14 


17.2 


.39 


400 


220 


1.30 


12.4 






15 


18.3 


• 35 


375 


220 


U15 


12. 2 






16 


19.4 


.31 


250 


210 


1.00 


12.1 






17 


18.3 


.33 


310 


200 


0.30 


12.1 






13 


13.3 


.45 


280 


170 


0.40 


12.0 






19 


17.2 


.52 


250 


140 


0.10 


12.0 






20 


16.1 


.59 


230 


190 


0.0 


11.9 





49 









TABLE : (Continued) 


. . , , 






Date 


EST 


(1) 
T airt°C> 


(2) 

Humidity 


(3) 
Wind 

Speed 


(4) 
Wind 

Direction 


(5) 
Solar 

Radiation 


(6) 

surf 


i 


21 


15,5 


.66 


210 


250 


0.0 


11.9 




22 


15.0 


• 72 


200 


310 


0.0 


11.3 




23 


15.0 


.70 


210 


340 


0.0 


11.8 




24 


15.0 


.68 


230 


10 


0.0 


11.8 


1/31 


1 


14.4 


.66 


250 


40 


0.0 


11.7 




o 


14.4 


.66 


250 


50 


0.0 


11.7 




3 


13. a 


.65 


270 


60 


0.0 


11.7 




4 


13.3 


.65 


290 


70 


0.0- 


11,7 




5 


12.7 


• ?7 


320 


90 


0.0 


11.7 




6 


12.2 


.89 


360 


110 


0.0 


11.7 




7 


11.6 


1.00 


400 


140 


0.0 


11.8 


1/31 


8 


12.7 


.97 


420 


130 


0.05 


11.3 




9 


14.4 


.95 


500 


220 


0.15 


11.3 




10 


16.1 


.93 


600 


250 


0.25 


11.9 




11 


17.2 


.36 


700 


230 


0.35 


12.0 




12 


16.6 


.79 


800 


310 


0.45 


12.1 




13 


16.1 


.72 


800 


330 


0.55 


10 




14 


15.0 


.74 


850 


330 


0.50 


12.1 




15 


13.8 


.76 


900 


330 


0.50 


11.9 




16 


12.7 


.77 


900 


330 


0.40 


11.7 




17 


11.6 


. 77 


7G0 


330 


0.30 


11.5 




IB 


10,5 


.77 


600 


330 


0.15 


11.3 




IS 


9.4 


.77 


600 


330 


0.05 


11.0 
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TABLE 2 [Continued} 






Date 


EST 


(1) 
T air(°C) 


(2) 

Humid itj 


(3) 

Wind 
Speed 


(4) 

Wind 
Direction 


(5) 
Solar 

Radiation 


(6) 

T 
surf 




20 


8.8 


.79 


600 


330 


0.0 


10.7 




21 


8.3 


.81 


600 


330 


0.0 


10.4 




22 


8.3 


.83 


550 


340 


0.0 


10,0 




23 


7.7 


.80 


500 


340 


0.0 


10.0 




24 


7.7 


.77 


500 


340 


0.0 


10.0 


2/1 


1 


7.2 


.74 


500 


340 


0.0 


10.0 




2 


6,7 


.74 


550 


340 


0.0~ 


10.0 




3 


6.1 


.73 


550 


350 


0.0 


9.9 




4 


5.5 


.73 


550 


350 


0.0 


9.9 




5 


4.4 


.73 


500 


350 


0.0 


9.8 




6 


3.8 


.73 


500 


360 


0.0 


9.8 




7 


3.3 


.73 


500 


360 


0.0 


9.7 




8 


3.8 


.68 


550 


360 


0.05 


9.8 




9 


5.0 


.58 


600 


360 


0.15 


9.9 




10 


6.1 


.51 


650 


360 


0.25 


10.0 




11 


7.2 


.45 


650 


350 


0.40 


10.1 




12 


8.3 


.39 


650 


350 


0.55 


10.2 




13 


9.4 


.33 


600 


350 


0.70 


10.2 




14 


10.0 


.30 


600 


340 


0.85 


10.1 




15 


10.5 


.32 


600 


340 


0.80 


10.1 




16 


11.1 


.34 


600 


330 


0.70 


10.0 




17 


10.0 


.40 


600 


330 


0.55 


10.0 




18 


8.3 


.46 


600 


330 


0.40 


10.0 
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Figure 1. Anclote Anchorage location in the state of Florida 
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Figure 2. Map of Anclote Anchorage 
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Fioure 3. Definition sketch of cr-coordinate 
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Figure 4. Grid arrangement in the horizontal proiection 

• (full-grid) depth (h) point; — u velocity point; 

! v velocity point; + (half-grid) w,a,p f T,o and. n point 
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(1+1.1+1) 




Mi.j) 3 



h(i+i.j+l) 



Figure 5. Four ceil* in a vertical column with velocities 
shown at definition point, and scalar variables 
at the center of cell M 
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Fiqur« 6. Notations (left) and variables (right) used in calculation 
of T (*), n m, u (** and v (*)* within f i ,11 -block 
(cross-hatched) 
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Figure 7. Grid work for the Anclote Anchorage 







Figure 8. Location of stations for in-situ measurement, June 1978 
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Figure 9. Velocity from In-situ measurement at 1710-1903, June 19, 1978 
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Figure 10. Velocity from in-situ measurement at 0648-0812, J 



une 20, 1978 
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Figure 11. Velocity from in-situ measurement at 1125-1245, J 




une 20, 1978 
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Figure 12. Velocity from in-sltu measurement at 1*50-1605, June 20, 1978 
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Figure 13. Daytime flight lines on June 19 and 20, 1978 
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Figure 14. Surface temperature in degree C from in-sltu measurement at 
1710-1903, June 19, 1978 * 







Fioure 15. Surface temperature in degree C from in-situ measurement at 
0648-0812, June 20, 1 ?78 







Figure 16. Surface temperature in degree C from in-situ measurement at 
1125-1215, June 20, 1978 
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Figure 18. Location of stations for In-sltu measurement. January 1979 
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Figure 19. Velocity from in-sltu measurement at 1020-1340, January 30, 1979 




Fioure 20. Velocity from in-situ measurement at 1 WO- 1800, January 30, 1 
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F!jure 21. Velocity from in-situ measurement at 1W0-1640, February 1, 1979 
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Fiqure 2 3* Surface temperature In degree C from In-sltu measurement at 
1M0-1800, January 30, ;979 
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Figure 21. Surface temperature In degree C from In-sltu measurement at 
1430-1640, February 1, 1979 







Figure 25. Surface temperature from In-situ measurement at 
1020-1340, January 30, 1979 
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Figure 26. Surface fl ten perature from In-sltu measurement at 
1110-7800, January 30, 1979 
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Figure 27. Surface temperature from )n-situ measurement at 
1430- 16*10, February 1 f 1979 
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Figure 28. Semidiurnal tide for June 19-20. 'l978 at south end of Anclote Key 
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Figure 29. Semidiurnal tide for January 30-February 1, 1980 at south end of Anclote Key 



TlflCtJUMC 20.1976H 1Q.S 

HI WO SfEE0(Cn/9eC)i 450.0 

UINO 0IRECTIONCO£O/N)i 110. 

AIR TCttFeRftTURC<OCG-C)i 26.6 

DISCHARGE TEf1f(0EG-C)i 29 .5 

0I3CH FtOMRRTE(CUn/3EC)i 62.7 

LENGTH 3CRLE(lCrt= X CM) i 41019. 

cLQCITT 3CflLE<Cn/3EU* 32.49 
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Figure 30. Surface velocity, Anclote Anchorage by mocielina 
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TIHEtJUME 20.1970)t 10. S 

UIMO Sf«0(CH/3tCli 450.0 

MI NO OlftECTIONCOeO/N)* 110. 

RIR TErtPCMTUfteCOeO-Clt 20 .0 

0I3CHA3OE TEflPCOEO-Cli 29 .3 

013CH FL0MHRTC(CUrt/3eCJ» €2.7 

LENGTH SCftLEUCfU X Cfl 1 i 41019. 

VELOCITY 3CRCE(C«/3EC)i 52.49 
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J* 4 
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Figure 31. UW velocity, Anclote Anchorage by modeling 
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TtriCC JUNE 20.1976)1 
H!N0 3PK0(CH/S€C»« 
HIND O!KECTION(OC0/M)t 
AIR TCnfCWTU«(0CO-CJi 
0I3CH*R0C TWP(0eO-C)i 
OUCH FtOHMTCCCUII/SeOi 
LCMOTM 3C«L2(1C*« X Cfl J • 
VCLQCITY 3CALECCt1/3CCli 



10.8 

480.0 

110. 

20.8 

29.5 

62.7 

41019. 

52.49 




Is 12 




Is 8 




U 4 



HIGH HOC 



Fiqure 32. VW velocity, Anclote Anchorage by modeling 
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TlnCCJUWe €0.1970)1 

mino sfceoccn/seoi 

MI NO OI«CTIONC0CO/Nli 
AIR TC«P€IWTURC(Oea-C)i 

ouchmoc TMfcoca-oi 

0I3CM FL0MWTCCCUn/3eC)i 

lcmotm scale ncn« X CflJi 

VELOCITY 3C«CC(Cf!/3£C>i 



11.3 

420.0 

no. 

29.1 
29.0 

02.7 

41019. 

32.49 
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Figure 33* Temperature from IR 
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TlfieCJUN* to. 1979 1 » 
UIHO 3fCCD(Cn/3CCJi 
MtMO 0l*eCTttMC0eO/N)i 

AtK Ttcnpcitnmitecoec^oi 

0I3CHWOC TEW (OCO-C) i 
0I3CM FU0H^ftTC(CUn/3CC)i 

lcmoth scflLencn* x cnn 

VELOCITY SCALE (C.V32C)i 
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Figure 3*t. Surface temperature, Anclote Anchorage by modeling 
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TIRECJUME 20.1978) t 14.8 

kino spccotcn/seoi 403.0 

H1N0 0IRECTI8M(0EO/N)i 110. 

AIR TEJirEBBTUREtOEO-O* . 80.8 

DISCHftftOE TEHfCOEO-CJi 80 .0 

OISCM EL8HRnTECCUn/8EC)i 62.7 

LENOTH SCRLEClCh* X Ct1)< 41019. 

VELOCITT 8Cm.ECCn/3ECli 82.49 
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Figure 3S. Surface velocity, Anclote Anchorage by modelinq 
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TlflElJUNC 20.1970)1 . 14.5 

H1N0 3KE0CCIt/3K)i. 400,0 

WINO OtMCTION(OeO/M)i 110. 

Rift TEflPEftftTUftEtOEOrOi 30.6 

0I3CHRR0E TCflftOeO-Oi 30. 

0I3CH FtOMftftTC(CUn/3CC)i 6E-7 

LENGTH 3CftLC(lCf1^X.Cri)t 41019. 

V6LQCITY 3CfltC(Cn/3CC)t 5E.49 
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J» 4 



E98 TIOC 



Figure 36. UW velocity, Anclote Anchorage by modeling 
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TMEUUNE 20.1970)1 14.9 

hi ho srccorcn/aecji 400.0 

UINO OIRECTiaN(0CO/N)i 110. 

RIR TEHPEMTUHEIOEO-Cli 30.6 

013CHRRCE TEHftOEO-Cli 30.0 

OtSCH FLOWARTECrUfl/SECJi 62. 7 

CEHOTM 3Cm.EC101s X Cfll t 41019. 

VELOCITY SCflLECCtt/SEUt 52.49 
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1= 8 




U 4 
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Figure 37. VW velocity, Anclote Anchorage by modeling 
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TitlECJUNE 2Q.1976)i 

uino spccotcn/secji 

MINO OIRECTIQN(0EO/N)t 
RIR TEnfERRTURE(0EO*C)i 
013CHM0E TEHPtOEO-CJi 
0I3CH FLQMftRTE(CUn/9EC)t 
LENGTH 3CALE(lCf1« X CHI • 
VELOCITY 3CALE(Cn/3EC)i 



15.0 
400 .0 
110* 
30.6 
29 .6 
62.7 
41019. 
52.49 




EB6 YIOE 



Figure 38. Temperature from IR 
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T met JUNE 20.1978)1 
HINO 3fCC0(CH/3CC)i 
UINO OIRECTIONCOEO/Hli 
RIR TEHPERATURECQEO-C)i 
(DISCHARGE TEflPtOEO-CJi 
0I3CH FtOMRATEtCUn/3ECJi 
LENGTH 3CALEtlCI1x X CMi 
VELOCITY 3CALE<Clt/3EC)t 



14.3 

400. Q 

no. 

30.6 
30 .0 
62.7 
41019. 
32.49 




DEVIATION PROn IR TEHP i 
EB8 TIDE 



0.361 



Figure 39. Surface temperature, Anclote Anchorage by modeling 
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TIMECJUNE 20.1978)1 

mino sr«o(cn/3coi 

MIHO OIRCCTIQN(OeO/N)i 
AIR TEHrER«TUREtOEO-C)i 
0I3CHAROE TEHPtOEO-Oi 
0I3CH PL0MRRTECCUI1/3EC). 
LENGTH 3CRIEUD1* X Cflli 
VELOCITY 3CRLECCrt/3EC)i 



17.5 

310.Q 

110* 

29.* 

30.2 

62.7 

41019. 

32.49 
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Figure 40. Surface velocity, Anclote Anchorage by Modeling 
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TttlECJUNE tO. 1976)i 

nino 3reco(cn/9ec)i 

NINO OIRECTIQN(QeO/N)» 
RIR TEnrEKATURECOEO-Cli 
DISCHARGE TCHf(0eO-C)i 
0I3CH FU5NKRTECCUH/3EC)i 
LENGTH 3CflUE(lCn» X CflJi 
VELOCITY 3C«LE(Ctt/3EC)i 



17*3 

310.0 

110. 

E9.4 

30. t 

62.7 

41019. 

32.49 





j» a 



J« 4 
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Figure 41. UW velocity, Anclote Anchorage by modelinq 
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TlftCMUNC 20.1979)1 17. S 

MINO 3fCE0(Cn/aCC)t 310.0 

MIMO 0!RCCTIQN(0EO/N)i 110. 

AIR TCHPM«TU«(0CO-Cn 29.4 

0I3CHAK0C TErtrtOCO-C)i 30.2 

013CH FL5HWTC(CU«/3CC)i 62.7 

LENGTH SCALSilCh. X Cm i 41019. 

VELOCITY 3CALE(Cfl/3EC)i 52.49 
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Ir 3 
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Figure 42. VW velocity/ Anclote Anchorage by modeling 
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TinECJUNC 20.1973H 
M1N0 3PECD(Cn/3€CJi 
NINO OI*tCTIQM(OeO/N)i 
AIK TCHPMRTURE(DCO-C)i 
013CHRH0E TCnP(OCO-C)i 
QI3CH PL0MR«TC(run/3CCJi 
LENGTH 3CflLEUCrt» % CfUi 
VELOCITY 3CRLC(Cn/3eC)i 



17.3 

310*0 

U0« 

29*4 

30.2 

62.7 

41019* 

32.49 




LQM TIDE 



Figure 43. Temperature from IR 
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TirtgfJUNE 20.1978)1 17,5 

mino spcEOtcn/seot 310.0 

MIND 0I*ECTI0N{0EO/N)i 110. 

AIR TEnPERRTUREfOEO-Cli 29.4 

0I3CHRROE TEttPtOEO-CJi 30. 2 

OUCH PL0WRRTetC\M/3EC>t 62.7 

LENGTH 3CALEilCria X CrtJi 41019. 

VELOCITY SCRL£{Cn/3ECJt 52.49 




OEVIflTtON FRQtt IR TEttr* 
L3H TIDE 



0.539 



Figure W. Surface temperature, Anclote Anchorage by modeling 
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IflCCJUNC 2Q.1976H 20*5 

WIND 3PCE0(Cn/3CC)t 225*0 

MINO OIRECTtONtQEO/Hli 90. 

IK TEHPERRTURECOEO-Oi ' 28. 

DI3CHRR0E TCrif(0EO-C)i 30.1 

0I3CM FLOHRRTEfCUh/SECli 62.7 

ENOTH 3CfllEtlCf1* X Cfl)t 41019. 

EL3CITY 3CRLC(Cn/3CC)i 52.49 
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Figure 45. Surface velocity, Anclote Anchorage by modeling 
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TIHEU'UNE 20.1976) i 

mino arccotcti/aeoi 

M1N0 0!ftECT!6N(CEO/N)i 

Htn Tcnrc«RTURC(oco-oi 

0I3CHMOE TEttftOEO-CH 
0I3CM *l5MRflTE(CUf1/3EC)» 
LENGTH SCALE ClCIs X Cfl 5 i 
VELOCITY 3CflLE(Cn/3EC)i 



20.5 

90. 

25*0 

30.1 

62.7 

41019. 

52.49 
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Figure 46. UW velocity, Anclote Anchorage by modeling 
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TtnCCJUNC 20*1978)1 

mino 9rceo(cn/aec)« 

MINQ OtKCCTI8H(020/>Mi 
AIR TCttfCR«TUKCC0CO-C)i 

dischrhoc TCHrtoca-oi 
o*3cm rLQMR^Teccun/aeoi 

LCNOTM SCAteUCn* X CfUi 
VELOCITY 3C«LC(Cn/3CC)t 



20*3 

225.0 

90* 

28*0 

30*1 

62.7 

41019, 

32*49 




ir 12 




U 8 




U 4 



RJJO HOC 



Figure 47* VW velocity, Anclote Anchorage by modeling 
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TIME1 JUNE 20.1978)1 
HINO SPEED* Cf1/3EC)i 
HINO 0IRECTIQN(0EO/M)i 
RIR TEf1PE*STURE<0EC-C)i 
DISCHARGE TEttPCOEC-Cji 
0I3CH PL3MRflTE(CUf1/3EC)i 
LENGTH SCfllEClCtls X Cf1)« 
VELOCITY SCflLECCf1/3ECJ* 



20.0 

225.0 

90. 

28.0 

29.3 

62.7 

41019. 

52.49 




FL30 TIOE 



Figure 48. Temperature from !R 



99 



TIMCC JUNC 20.1970)1 
NINO 3rCC0(Cfi/3CC)i 
NINO 0XRECTIQN(0EO/N)i 
fl!R TEMPCR«TURC(06O-C)i 
0I3CHAKOC TCriP(0CO-C)» 
0I3CH FLQNRRTE(CUrt/3CC)i 
LCNOTH 3CflL£(lCfl3 X CH)i 
VELOCITY 3CflL£tCn/3EC)» 



20.5 
223.0 

90. 

20.0 

30.1 

$2.7 

41019. 

52.49 




DEVIATION PRQft IK TErtf i 
FLOO TiOC 



0.361 



Figure 49. Surface temperature, Ancfote Anchorage by modeling 
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TlhEf JRNURRY 30*1970 J : 

MtNO 3PEE0(Cf1/3ECJi 
WtNQ OiRECTlONfGEO/Nj: 
AIR TEttPERfHUREfOEG-CJ : 
OlSCHfiRCE TENPfOEO-CJ: 
OtSCH FLOHRflTtrCUM/SECJ i 
LENOTH 3CflLE(lCri* X Cflh 
VELOCITY 3CALE(Cri/SEC)» 



11.0 

400.0 

320. 

12.5 

15 c 6 

62.7 

41019, 

52-49 
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Figure 50. Surface velocity, Anclote Anchorage by modeling 
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UMEf JANUARY 30-1979)1 

wino 3pecorcfi/3to» 

HINO 01RCCT10NCCE0/N)i 
AIR UftPERATUREfOEO-Cii 
DISCHARGE TfcflPCOEO-CJi 
015CH PLOMRATEf CUM/SEC n 
LENOTH SCflLEftCris X CM ) « 
VELOCITY 3CflUfCfl/3EC)> 



n.o 

400*0 

320. 

12.5 

IS. 6 

G2.7 

41019. 

52.49 




j* a 




J* 4 



KQO 710E 



Figure 51. UW velocity, Anclote Anchorage bv modelina 
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T r t £(JflNUflRr 30. 1979 i * 
WIND SPEEOCCN/SECJi 
NINO OIRECTlONfOEO/NJi 
fllft TEftPERfiTUREfOEC-CJi 
DISCHARGE TFMP(0EO-CJ» 
OiSCH FL0WRRUfCUff/5EC)» 
LENGTH SCPLEflCMr X Cf^* 
VELQCUr 5C.U£rCM/5£CJi 



400*0 
320. * 
12.5 
IS. 6 

62.7 

41013, 

b<:.4i3 




1= 12 




i= a 




r 4 



FLOO TIOE 



Ficure 52. VW velocity, Anciote Anchoraae by moricling 
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TlttEf JflNURRY 30.19791* 
HINO 3rEE0CCrt/3tCJ> 
NINO OIRECTIONCOEO/N). 
«IR TEriPERATUREfOEO-C)> 
013CHR*CE TEflPtOEO-CJi 
013CH PL0WRRTEfCUtt/3ECn 
LENGTH SCRLEflCrir X CflJ • 
VELOCITY 3CALE<Crt/3ECJi 



n.s 

400.0 

320. 

12.7 

15.5 

62.7 

41019. 

52.49 




FLOO T10E 



Figure 53* Temperature from 1R 
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UftEfJPNURRf 30,1979)1 
MNO SPEEO(C«/SEC)i 
WINO OIRECTIQN(OEO/N)i 
PIR TErtPERPTUREfDEO-C)* 
DISCHARGE TEP1P(GE0-C)» 
OiSCH FLOHRATEfCUfl/SECii 
LENGTH SCPlEflCttr A CM ) i 
VELOCITY SCflLEf CM/SEC J t 



U.Q 
400 .0 
J20. 
U.5 
IS. 6 
G2-7 
41019, 
52.49 




DEVIATION FROft IR TEflPt 
5100 TIOE 



0*742 



Figure 54. Surface temperature, Anciote Anchorage by modeling 
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tlflEf JANUARY 30.1979) « 

mino 9PEE0f cm/sec )« 

HINO OIRECTlONrOCO/NH 
AIR TEf1PERA7UREfOEO-C)i 
0I9CHARCE TEflPfOEO-CJi 
019CH FLQMRA1ErCUf1/9ECJ> 
LtNCTH 9CAlEflCf1r X Cfl ) t 
VELOCITY 9CALEftf1/9ECJi 



16-0 

400.0 

360. 

H.O 

1G.J 

62.7 

41019. 

52.49 



\ V 



t i 



' * ' ' S S S S / / ' ' ' 

' / / / / S ' ' " " 
1 I / / / '■ ' ' 

W / / / ' ' 

I / / / ' ' ' 

n ' / / / / ' ' ' 

'lilt'" 

I / / / ' ' ' 



tflB TlOfc 



Figure 55. Surface velocity, Anclote Anchorage by modeling 
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Tintr JPNOPRY J0.t979Ji 
H1NQ 3PEE0fCN/SEC)i 
U1NO 01PECT!QN(OEO/NJi 
AIR TENPEKfmjREfOEO-CJi 
0I5CHPRCE TtHP(OEO-CJi 
DiSCri FLOMRRtEfCUH/SECJi 
LENGTH SCPLcf tCNr X CM J 1 
VELOCnr 3CaiEfCM/5EC]« 



ICO 

400*0 
360* 
14.0 
16.3 

62.7 

41013. 
5^.43 



J* 12 




Jr d 




Jr 4 



£38 net 



Figure 56. UW velocity, Anclote Anchorage by modelinq 
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UrtEfJflhUPKf 30.1079J » 
M1N0 SPEEOf CM/5CCJ i 
MINO QlRCCTlONf OCO/N)» 

aia TtfireKRiuReroto-cji 

DISCHARGE TMPfOtO-CM 
013CH FlOM««UfCUH/«CJ« 
LENGTH SCPLtflCrU X Crt ) i 

veiocm acflitfcn/atcn 



ICQ 

400.0 

560* 

14.0 

16.3 

41013. 
52.49 




U 12 




lr 8 




Ir 4 



Eae Tioe 



Figure 57. VW velocity, Anclote Anchorage by mocielinq 
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7IMEURNUPRY 30.1973 J > 
U1N0 SPEE0(CM/3EC)> 
U1N0 01RECT!QN(0eO/N)i 
FIK UttPERflTURE(QEO-C)t 
013CHPRGE TEHPCOEO-CM 
015CH FiaWRAUfClifl/SEOi 
LENGTH SCflLEflCflr X CH ) » 
VELOCITY 3CALE(Cf1/3ECJ> 



n.o 

400*0 
360. 

ie.0 

16.3 
62.7 
41019. 
52.49 
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Figure 58. Temperature from 1R 
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TIfltf JANUPAY JQ.t370)< 
WiHO 3f«0(Cn/UC)» 
UiNO OtRCCTlOH(OeO/N)i 

rih unrert«iuRtroeo-c>i 
oiscmauge unrcoeo-u* 

015CH PlQMRAU(CUH/3EC). 

itNOTH scaur ten* x cn» i 

VfcLOCPY 3CAUftn/3tCJi 



ig-o 

400*0 

J6Q. 

14.0 

tea 

6* .7 
41010- 
Si* 40 




CtVlftUOH MOri IK UftFt 

tsa not 



0*G50 



Figure 59. Surface temperature, Anciote Anchorage by modeling 
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TIflEf 1ST FEB. .137*3 i 
WIND SPEEOrCtt/SCC J i 
UINO OIRECTtQN(0EO/M>t 
RIR TEnrERATURE(OEG-C)i 
DISCHARGE TeMPtOCC-C]! 
013CH FLflWRRT£(CUn/3EC! i 
LENGTH SC«LE( ICHr X CM ) i 
VELOCITY 3CALEtCM/SECl i 



16.0 

600. 

330^0 

11 .1 

M.O 

62,7 

41013, 

52.49 




HIGH TiQ£ 



Figure 60* Temperature from 1R 
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